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PREFACE 


This  technical  report  comprises  publications  on  the  reflection  of 
ultrasonic  waves  at  a liquid-solid  interface  made  between  January  1975 
and  the  present.  During  this  period  a number  of  publications  have  been 
made  on  other  subjects;  however,  we  prefer  to  concentrate  on  the 
reflection  problem  in  this  technical  report  in  order  to  present  a 
unified  picture  of  the  development  of  this  subject.  A later  technical 
report  will  cover  other  subjects. 

The  report  begins  with  the  English  language  version  of  a paper 
whose  Russian-language  counterpart  was  included  in  Technical  Report 
No.  11,  It  may  be  worthwhile  to  point  out  that  our  Russian  colleagues 
were  impressed  with  the  paper  enough  to  make  it  the  lead  article  in 
Akusticheskii  Zhurnal  for  the  year  1975. 

In  this  paper  we  were  able  to  get  agreement  between  theory  and 
experiment  at  the  Rayleigh  angle  sufficient  to  permit  us  to  plot  them 
on  the  same  graph.  (Up  until  this  time  theory  and  experiment  were 
agreeing  only  qualitatively  so  that  such  a comparison  was  not  possible.) 
The  agreement  was  made  possible  by  advances  both  in  experiment  and  in 
theory.  The  experimental  advance  was  the  use  of  a transducer  which 
produced  a truly  Gaussian  amplitude  distribution  in  the  incident  beam 
(see  Technical  Report  No.  8).  The  theoretical  advance  was  the  retention 
of  higher-order  terms  in  the  power  series  expansion  of  the  phase  shift 
upon  reflection  than  had  been  used  previously  by  Brekhovskikh.  This 
led  to  a complicated  expression  for  the  reflected  beam,  a power  series 
expansion  in  Fresnel  integrals,  which  had  to  be  evaluated  by  computer. 


ii 


iii 

A mathematical  improvement  resulted  from  interaction  between  the 
author  and  Werner  G.  Neubauer  and  Larry  Flax  of  NRL.  Rather  than  a 
series  expansion,  a closed-form  solution  in  terms  of  error  functions 
was  introduced.  This  is  contained  in  the  second  paper. 

In  the  third  paper  we  show  that  a considerable  improvement  of  the 
agreement  between  theory  and  experiment  results  from  the  use  of  the 
theory  of  Bertoni  and  Tamir,  with  modification  to  account  for  the 
propagation  distance  of  the  incident  beam.  The  data  in  this  paper  also 
are  very  carefully  taken. 

The  fourth  paper,  given  at  the  XV  International  Conference  on 
Acoustics-Ultrasound  in  Prague,  is  a summary  paper  describing  the 
application  of  the  theory  of  Bertoni  and  Tamir  to  the  results  on 
stainless  steel.  In  this  paper  we  also  indicate  that  a corrugated 
interface  makes  possible  coupling  to  a surface  wave  which  propagates 
in  the  negative-x  direction  along  the  interface.  The  fifth  paper  shows 
that  this  negatively  directed  surface  wave  causes  the  reflected  beam  to 
be  displaced  in  the  opposite  direction  from  that  observed  at  the 
Rayleigh  angle. 

The  sixth  paper  shows  not  only  that  the  displacement  observed  at 
liquid-solid  interfaces  can  be  observed  at  water-sediment  interfaces, 
but  also  that  the  backward  displacement  phenomenon  might  exist  because 
of  the  periodicity  resulting  from  the  granularity  of  the  sediment. 
Furthermore,  the  water-sediment  interface  corresponds  to  a water-solid 
interface  for  which  > V > Vg,  a situation  which  previously  had  not 
been  observed  to  give  a displacement  of  the  reflected  beam.  However, 
such  a displacement  does  exist  for  this  class  of  liquid-solid  reflectors. 
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as  we  show  in  the  seventh  paper.  The  displacement  is  observed  near  the 
critical  angle  for  water-plexiglass,  not  the  Rayleigh  angle,  however, 
for  there  is  no  Rayleigh  angle. 

The  eighth  paper  is  an  application  of  the  surface  wave  phenomena  to 
the  detection  of  subsurface  flaws  in  metals.  The  surface  waves  generated 
at  the  Rayleigh  angle  are  described  in  the  theory  of  Bertoni  and  Tamir 
as  "leaky  waves"  because  the  energy  "leaks"  from  them  back  into  the 
liquid  to  become  the  reflected  beam.  The  ninth  paper  shows  how  these 
leaky  waves  can  be  used  to  detect  subsurface  flaws  in  solids.  The  final 
paper  begins  a study  of  the  effect  of  the  introduction  of  a layer  on  the 
solid. 


M.  A.  Breazeale 
February  1978 
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Paper  No.  1 


Energy  Redistribution  of  a Gaussian  Ultrasonic  Beam  Reflected  from 
Liquid-Solid  Interface  (M.  A.  Breazeale,  Laszlo  Adler,  and 
James  H.  Smith),  Akust.  Zh.,  21,  1-10  (1975)  [Sov.  Phys.  Acoust. 
21^,  1-6  (1975)]. 


[2] 


Energy  redistribution  of  a Gaussian  ultrasonic  beam  reflected  from  liquid— solid  interface 

M.  A.  Breazeale,  Laszio  Adler,  and  James  H.  Smith 
Department  of  Phyilci,  Untvenlty  of  Tenneiiee.  Knoxville,  Tenncuee 
(Submitted  February  4,  1974) 

Akust.  Zb.,  ^ 1-10  (Januaty-Febnury  1975) 

The  problem  of  the  reflection  of  an  ultrasonic  beam  from  a liquid— solid  Interface  Is  considered.  The  concept 
of  'beam  displacement*  Is  reexamined  by  studying  the  Iniernal  structure  of  the  reflected  beam.  This  study  Is 
carried  out  by  meant  of  a transducer  which  produces  a Gaunlan  distribution  of  energy  across  the  beam  width. 

Experimental  results  show  that  the  energy  of  the  reflected  beam  It  redistributed  Into  two  or  more  beams  for 
angles  neat  the  angle  of  optimum  generation  of  surface  waves  on  the  Interface.  Reasonable  agreement  is  ob- 
tained between  theory  and  experiment  for  water— aluminum  and  water— brass  Interfaces. 


1.  INTRODUCTION 

Consider  an  ultrasonic  beam  incident  on  a liquid—  solid 
interface.  For  a wide  range  of  angles  of  incidence  this 
problem  can  be  mathematically  described  by  solving  the 
boundary- value  problem  under  the  assumption  of  energy 
conservation.  At  the  critical  angles,  however,  experi- 
mental verification  of  the  calculated  results  is  compli- 
cated by  the  fact  that  the  reflected  beam  appears  to  be 
"displaced."  This  effect  is  especially  pronounced  at  the 
angle  for  excitation  of  surface  waves  along  the  interface, 
as  has  been  demonstrated  by  Schoch.*  The  "displacement" 
of  the  reflected  beam  is  attributed  to  the  fact  that  a beam 
does  not  behave  as  an  infinitely  extended  plane  wave,  as 
is  tacitly  assumed  in  the  simplified  theory. 

The  situation  is  shown  in  Fig.  1 . An  ultrasonic  beam 
is  incident  upon  an  interface  at  an  angle  9g  (approximately 
the  angle  at  which  surface  waves  are  excited  on  the  in- 
terface). The  beam,  assumed  to  be  Infinitely  extended  in  a 
direction  normal  to  the  figure,  is  displaced  a distance  A 
upon  reflection.  In  developing  a theory  to  describe  this 
effect,  one  includes  the  mutual  phase  relation  of  partial 
waves  in  the  reflected  beam  by  expanding  the  expression 
for  the  phase  shift  upon  reflection  4>(p)  into  a power  series 
of  the  form 

d>(y)  (p-a)  + 'lt<J>" (p—o)' ....  (1) 

where 


y— A sin  0,  a— fc  sin  6i,  At— 2n/X, 


X = the  ultrasonic  wavelength  in  the  liquid,  Sg  is  the  angle 
of  incidence,  and  6 is  the  angle  of  incidence  of  a wavelet 
in  the  Incident  beam.  By  carrying  out  this  analysis  in  the 
first  approximation,  one  finds  that  the  displacement  is 


In  this  approximation,  the  beam  profile  of  the  reflected 
beam  (l.e.,  the  amplitude  distribution  of  the  reOected 
beam  measured  parallel  to  the  Interface)  is  the  same  as 
the  beam  profile  of  the  Incident  beam.  Such  profiles  are 
sketched  on  Fig.  1 as  Fig.  la  and  Fig^ure  lb. 

In  making  a comparison  between  theory  and  experi- 
ment, it  has  been  found  that  there  is  nominal  agreement 
between  the  measured  "displacement*  and  that  calculated 
from  relation  (1);  however,  the  e}q>eriment  is  compli- 
cated by  the  presence  of  unexpected  reflected  beams.  For 
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this  reason,  it  is  necessary  to  carry  out  the  analysis  in 
the  second  approximation,  to  include  the  second-order 
terms  in  relation  (1) . In  this  higher  approximation,  the 
"beam  displacement"  loses  Its  meaning.  The  profile  of 
the  reflected  beam  is  not  the  same  as  that  of  the  incident 
beam,  l.e.,  there  is  an  "energy  redistribution"  in  the  re- 
flected beam.  (Figure  lb  would  no  longer  be  the  same 
as  Figure  la.) 

Brekhovskikh  * has  developed  a detailed  theory  for 
the  reflection  of  a bounded  ultrasonic  beam  at  a liquid- 
solid  Interface.  From  this  theory  he  obtained  an  expres- 
sion for  the  amplitude  of  the  reflected  beam  and  for  the 
■displacement"  of  the  reflected  beam.  Brekhovskikh 
also  pointed  out  that  the  structure  of  the  reflected  beam 
Is  in  general  different  from  that  of  the  incident  beam,  l.e., 
that  there  is  an  "energy  redistribution."  He  obtained  an 
expression  for  the  amplitude  distribution  in  the  reflected 
beam  as  a function  of  the  distance  along  the  interface,  and 
solved  a special  case  under  the  assumption  that  the  in- 
cident beam  is  accurately  described  by  a step  function. 

In  a comparison  between  theory  and  experiment,  we 
have  encountered  the  fact  that  a step-function  distribution 
In  the  incident  ultrasonic  beam  leads  to  diffraction  effects 
which  confuse  the  comparison.  For  this  reason,  we  have 
developed  a transducer  which  produces  an  amplitude  dis- 
tribution in  the  incident  beam  which  is  accurately  described 
by  a Gaussian  function.^  The  Gaussian  function  amplitude 
distribution  produces  a single  well-defined  beam  without 
diffraction  lobes.  The  purpose  of  the  present  work  is  to 
show  how  closely  the  e3q>erimental  results  obtained  with  a 
Gaussian  distribution  of  amplitudes  in  the  incident  ultra- 
sonic beam  agree  with  those  calculated  from  the  approx- 
imate extension  of  the  Brekhovskikh  theory. 

2.  ENERGY  REDISTRIBUTION  OF  A 
REFLECTED  GAUSSIAN  BEAM: 

THEORY 

In  order  to  represent  the  beam  as  a superposition  of 
plane  waves,  Brekhovskikh  writes  the  field  in  the  plane 
z = I (at  the  interface)  as  a Fourier  integral.  By  Intro- 
ducing the  amplitude  reflection  coefficient  V(9j),  which  is 
a function  of  the  angle  of  incidence  9g,  the  general  form 
of  tlie  solution  to  the  wave  equation  which  describes  the 
reflected  beam  is  written  as 

♦■ten  W— -^^*'(0.)e*P{iIox+/(*^-a*)  IJ) 

u*)jF(i+<«>'+yn<f>''u)du,  (3) 

Copyrighi  01975  Ameilcxn  Inrtitule  of  Phydc* 


1 


1 


[3] 


rig.  1.  Fint-order  repceaentatlon  of  the  reflection  of  in  ultriionic  wive  it 
I llquld'-solid  interfice.  i)  Incident  beam  cros  section:  b)  reflected  beim 
cron  section. 
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+r,+r.+7-.+r,+r,)+(r.+7',+7’„+r„+r,.+r.,+r»)*i-  (8) 

where 

r,-C(B.)-C(u.);  (9) 

T,—2B'0"'  Juj'sin  ("Y"*’) 

-nu,*  sin  u,*  J -2  cos  o,"  J J d®) 


where  I is  the  normal  distance  between  the  transmitter 
and  the  liquid— solid  Interface,  u = (x  + *'  — 

Is  the  variable  of  integration  in  the  original  Fourier  in- 
tegral, and  F(x  + *'  + u)  describes  the  incident  beam 
Amctlonal  dependence  on  the  amplitude  distribution  along 
the  Interface.  (In  Eq.  (3)  the  higher  approximation,  in- 
clusion of  second-order  terms  in  Eqi  (1),  has  been  made.] 

For  the  present  purposes,  F(x+*'  + u)  is  to  be 
a Gaussian  function.  Such  a function  does  not  change  ap- 
preciably over  a distance  of  one  wavelength  along  the 
propagation  direction.  This  satisfies  the  requirement 
Imposed  by  the  manner  in  which  Eq.  (3)  was  derived. 

The  Gaussian  function  describing  the  incident  beam  can 
be  written 

P(x+0’+  tn^a)-A  eip[-i?(*+®'+ynd>'o)*J,  (4) 

where  A and  B are  experimentally  determined  constants 
which  depend  on  the  characteristics  of  the  transducer. 
The  amplitude  of  the  reflected  beam  then  may  be  obtained 
by  substituting  Eq.  4 into  Eq.  3: 


4’ten ^ {‘l«+ J) 

XA  u*  j]eipI-B(i+d>'+yn<P"tt)*]}du.  (5) 

A solution  to  Eq.  (5)  can  be  obtained  if  the  limits  on 
the  integrals  are  defined  as  follows: 


— u,  sin  o,*  J +5  (u.)  J ; (11) 

r,— d>"B*^2nu,*»in  J +-^u,cos^-^u,) 

^C(u,)— 2nu,’sin /— a,)  — ^o,cos(  — u,\  +— =-C(u,)l; 

12  \2  f 2 \ 2 ! 12  ■» 

(13) 

r. — i^(x+0-)B‘  [ y “>'»■“(  Y “•)  * +0®  ( Y “•) 


(14) 

r, d>"'B*(x-t-(P')*  j2.5u,sm|-^ii,j  -2nSM 

— 2.5u,  sin  ^-^u,  J -i-2n5(u,)J; 

(15) 

7-.-S(u.)-5(u.); 

(16) 

T, 2B'0"'  [-nu.-cos  (y“*') 

+nu,*cos  (y“>')“28in^YU,*JJ  , 

(17) 

Ui~  — ■ (a+z+<P');  (8) 

1n0“ 

and 

at—  + —^=(a~x-0'),  (7) 

Tnth 

where  a is  the  half-width  of  the  Gaussian  beam  at  a point 
where  the  energy  decreases  to  l/e  of  its  maximum  value. 
Then,  using  e"*!  “ cosy  — 1 sin  y to  express  the  first  ex- 
ponential under  the  integral,  and  expanding  the  second  ex- 
ponential into  a power  series,  one  obtains  a series  whose 
terms  are  integrable.*  Terms  up  to  fourth  order  are  kept 
in  the  following  expressions.  The  energy  distribution  for 
the  reflected  beam  is  then  obtained  by  multiplying  the 
integrated  form  of  Eq.  (5)  by  its  complex  conjugate: 


r..— 8B*yy<P"*(x+<J>')  [-II.C08(yu/)  +C(u.) 

+u,  cos  u,*  j —C  (u.)  J ; (18) 

r„— 4B’(P''(i+<t>')*  [-«>»(y“*’) 

r,,— J-2nu,*cos  Y “•*'“( T””*  ) 

^S(ui)— 2nu,*cot  (— u,\  +-^u,8ui  (— u.)  +-^5(u,)1: 

12  \2  ' 2 12  / 12  ■* 

(20) 
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14] 


r„— ^l.lu,*j>ia  (yUt)  +0.9co3^YU,) 


-l.lUi’sin  ^-^u,  J -0.9cos^-^u,  J J ; 

T,, 0”H‘(xiO')  j^-2.5u.cos  (-Y“«)  +2nC(ii,) 

+2.5u,  cos  Ui  j -2nC(u,)  J ; 


(21) 


(22) 


where  the  C(u)  and  S(u)  are  tlie  cosine  and  sine  Fresnel 
Intcijrals,  respectively.  The  other  parameters  have  been 
determined  experimentally,  and  Eq.  (8)  has  been  program- 
med for  a computer  calculation. 


3.  ENERGY  REDISTRIBUTION  OF  A 

REFLECTED  GAUSSIAN  BEAM: 

EXPERIMENT 

Apparatus 

The  experimental  apparatus  is  shown  in  Fig.  2.  The 
electronic  system  is  a standard  pulse  arrangement,  with 
tlic  signal  from  the  receiver  transducer  being  displayed 
on  an  oscilloscope.  Two  aspects  of  the  arrangement  are 
unique,  however:  the  goniometer  and  the  transducers. 

In  the  experimental  situation  both  the  driver  and  the  re- 
ceiver transducer  are  mounted  on  the  goniometer  and  this 
unit  is  placed  in  a tank  of  water. 

Goniometer  system.  A schematic  diagram  of 
the  goniometer  is  shown  in  Fig.  2.  The  system  is  de- 
signed so  that  the  face  of  the  specimen  (or  reflector)  forms 
a vertical  plane  at  the  eenter  of  a circle  formed  by  rota- 
tion of  the  tw'o  arms.  A vernier  scale  on  each  arm  permits 
the  angular  position  of  each  arm  to  be  read  to  two  minutes 
of  arc.  Both  arms  have  transducer  holder  stages  mounted 
on  riders  that  allow  the  transducers  to  be  moved  up  to 
50  cm  in  a radial  direction.  In  addition,  the  transducer 
on  the  receiver  arm  can  be  moved  perpendicular  to  the 
arm  (or  radius).  This  way  accurate  scanning  of  the  re- 
flected beam,  can  be  accomplished.  The  reflectors  have 
dimensions  of  10x5x2.5  cm  and  have  machined  surfaces. 

Transducer.  Details  of  the  behavior  of  an  ultra- 
sonic wave  reflected  from  an  Interface  depend  fundamen- 
tally on  the  energy  distribution  in  the  cross  section  of  the 
incident  ultrasonic  beam  (see  Eq.  (3)].  The  mathematical 
interpretation  of  the  results  Is  facilitated  if  the  energy 
distribution  in  the  incident  beam  cross  section  can  be  de- 
scribed in  a simple  relation  such  as  a Gaussian  function. 
The  experimental  measurements  are  also  more  straight- 
forward if  one  uses  a Gaussian  function  distribution  in 
the  ultrasonic  beam.  This  can  be  demonstrated  by  schlieren 
photography. 

Flgurc3a  is  a schlieren  photographof  a 4 MHz  ultra- 
sonic beam  from  a narrow  transducer  incident  at  an  angle 
of  35  deg  upon  a water— aluminum  interface.  This  bc.am 
has  a step  function  distribution  of  amplitudes  near  the 
transducer.  The  interpretation  problem  resulting  from 
the  diffraction  lobes  is  obvious.  Figure  3b  shows  an  ul- 
trasonic beam  of  approximately  the  same  width  and  fre- 
•Juency  with  a Gaussian  amplitude  dl.strlbution  across  its 


width.  This  beam,  produced  by  a new  tnmsduccr,’  l.s  in- 
cident on  the  same  interface  at  the  same  angle  of  inci- 
dence. In  addition  to  the  abst.ice  of  diffraction  lobes, 
there  is  less  divergence  In  the  Gaussian  ultrasonic  beam. 

The  Gaussian  distribution  was  obtained  by  using  a 
narrow  electrode  at  the  back  of  the  quartz  transducer  and 
grounding  the  entire  front  plane.  Fringing  of  the  elec- 
trical field  in  the  quartz  transducer  material  closely  ap- 
proximates a Gaussian  function  if  the  ratio  of  the  width 
W of  the  back  electrode  to  the  thickness  T of  the  trans- 
ducer material  Is  2 < W/T  < 4. 

The  ultrasonic  beam  profile  emitted  by  these  trans- 
ducers could  be  obtained  by  using  the  goniometer  system 
with  the  sample  removed.  All  mea.surements  wilh  the 
goniometer  system  were  made  at  a frequency  of  2 MHz. 

A plot  of  the  energy  distribution  obtained  by  using  the 
new  transducers  both  as  transmitter  and  receiver  is  given 
In  Fig.  4.  The  experimental  points  have  been  fitted  by  a 
Gaussian  function  of  the  form , A exp(—  Bx^),  where  A = 1 


Fig.  2.  Experimental  apparatus  (dia- 
gram of  goniometer.  1)  Specimen(re- 
flector):  2,  3)  rotating  arms;  4)  mov- 
able emitting  transducer;  5)  movable 
receiving  transducer. 


Fig.  3.  Schlieren  photographs  of  transducer  outputs,  a)  Step-functlop  out- 
put; b)  Gauulan  output. 
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Fig.  4.  Cron  seciTon  of  transducrr  output. 


Fig.  5.  ScUleren  photogtaphs  of  actual  energy  ledlitilbutlon  In  reflected 
Gaunfaa  ultrasonic  beam,  a)  t^<  Op;  b)  0^  » dp.  the  angle  for  excitation 
of  Rayleigh  surface  waves;  c)  Oj  > dp. 


and  B = O.OIG.  The  agreement  between  the  calculated 
curve  and  the  experimental  points  Indicates  that  the  ultra- 
sonic beam  produced  by  tlie  new  transducers  Is  accurately 
described  by  a Gaussian  function. 

Experimental  Procedure 

To  be  able  to  compare  experimental  results  on  the 
energy  distribution  in  the  reflected  beam  with  the  the- 
oretical predictions,  It  is  necessary  to  determine  both  di' 
and  , In  addition  to  determining  the  slope  of  the  mea- 
sured and  calculated  phase  shift  <1>  done  previously,*  one 
can  determine  directly  by  measuring  the  "displace- 
ment," and  compare  these  results  with  those  calculated 
from  the  theory  of  Brekhovskikh.  The  quantity  d’"  is  then 
the  slope  of  d-'  at  the  appropriate  angle  of  incidence. 

Measurement  of  d»'.  From  relation  (2),  the 
quantity  d>'  is  the  negative  of  the  "displacement."  This 
Is  measured  as  follows.  The  transmitter  arm  of  the  go- 
niometer is  set  at  the  desired  angle  of  incidence  and  the 
receiver  arm  is  rotated  until  a maximum  siipial  is  ob- 
tained on  the  oscilloscope.  The  receiver  transducer  Is 
then  shifted  laterally  (perpendicular  to  the  receiver  arm) 
In  both  directions.  An  increase  of  received  signal  in- 
dicates the  direction  of  the  "displacement."  A maximum 
In  the  received  signal  locates  the  position  of  the  displaced 
beam.  This  displacement  Is  then  measured  for  different 
angles  of  incidence.  A correction  is  made  for  the  fact 
that  the  motion  is  perpendicular  to  a radius  rather  than 
parallel  to  the  Interface,  as  assumed  in  the  theory. 

Measurement  of  the  energy  redistribu- 
tion ■ The  energy  redistribution  is  more  meaningfully 
studied  by  measuring  the  energy  distribution  in  the  re- 
flected beam  both  at  the  angle  at  which  surface  waves 
are  excited  on  the  interface  (maximum  displacement)  and 
at  angles  at  which  displaceme:<^  is  imperceptible.  The 
procedure  was  as  follows. 

By  means  of  the  goniometer,  the  incident  beam  was 
set  at  a previously  chosen  angle  of  Incidence.  The  re- 
ceiver arm  was  set  at  the  same  (reflection)  angle.  The 
receiver  transducer  was  then  displaced  perpendicular  to 
the  receiver  arm  in  increments  of  2 mm  in  both  direc- 
tions. The  receiver  transducer  output  was  monitored  by 
means  of  the  oscilloscope.  Considerable  care  was  neces- 
sary to  determine  the  exact  angle  of  Incidence  for  max- 
imum displacement. 


rig.  <.  CroM  tection  of  leflected  GiussUn  beimi  it  0^  < 6^(i)  and  0^  > 

egfb). 


4.  RESULTS  AND  DISCUSSION 

Careful  schlleren  photography  has  shown  us  that  even 
with  the  transducers  which  produce  a Gaussian  distribu- 
tion of  energy  in  the  incident  beam,  the  beam  Is  reflected 
as  multiple  beams  when  the  angle  of  Incidence  is  that  at 
which  surface  waves  are  excited  on  the  Interface  (max- 
imum displacement).  This  is  shown  In  Fig.  Sb  for  a 
water— brass  Interface.  Figs.SaandSc  are  schlleren 
photographs  for  angles  of  Incidence  respectively  smaller 
and  larger  than  this  angle. 

Data  from  scans  of  the  reflected  beams  with  the  go- 
niometer system  for  angles  corresponding  to  Figs . 5a  and 
So  are  shown  in  Fig.  C.  A brass  reflector  was  used.  An 
Incident  angle  of  10*  Is  well  below  the  longitudinal  critical 
angle  for  brass;  an  angle  of  54*  Is  well  above  the  angle 
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for  excitation  of  surface  waves.  A stralf^htforward  cal- 
culation from  the  simplified  theory  gives  an  energy  re- 
flection coefficient  of  801  at  10*  and  100%  at  54*.  These 
values  are  consistent  with  the  maximum  values  of  the 
curves  In  Figure  6. 

« 

The  displacement  (of  the  reflected  beam  in  a direc- 
tion parallel  to  tlie  Interface)  as  a function  of  the  angle  of 
Incidence  was  determined  for  a water— aluminum  inter- 
face and  a water- brass  Interface.  The  results  are  shown 
in  Figs. 7 and  8.  Experimental  data  are  indicated  by  points. 
The  curves  arc  the  results  of  computer  calculations  based 
on  the  theory  of  Brekhovskikh.  On  these  curves  the  val- 
ues of  the  critical  ai\gle  for  the  longitudinal  wave  In  the 
solid  (Oqi)  and  for  the  shear  wave  In  the  solid  (Ocs)t  as 
well  as  the  optimum  angle  for  excitation  of  surface  waves 
on  the  interface  (Op)  arc  indicated.  The  maximum  dis- 
placement of  the  pi'incipal  beam  occurs  at  (?p.  Also,  at 
fip  there  Is  an  energy  redistribution  as  shown  in  Fig.  5b. 

The  graphs  in  Figs . 7 and  8 are  plots  of  the  negative 
of  the  quantity  4>',  as  indicated  in  Eq.  (2).  That  value  of 
most  Importance  to  the  following  comparison  between  the- 
ory and  experiment  Is  the  maximum  value  occurring  at 
Op  = 30.5”  for  aluminum  and  flp  = 47.0°  for  brass.  The 
values  of  these  quantities  evaluated  from  these  data  and 
used  in  the  computer  calculation  of  Eq.  (8)  were  (“flp  = 
-14.5  mm  for  aluminum  and  (4>')r  =-16.5  mm  for  brass. 

The  quantity  4‘''  could  also  be  determined  from  these 
curves,  although  this  was  a much  less  precise  measure- 
ment. The  values  used  in  the  computer  calculation  of 
Eq.  (8)  were  (4>’’)p  = — 5 mm*  for  aluminum  and  (4<")p  - 
— 7 mm*  for  brass. 

Energy  Bedlstributlon  in  the 
Presence  of  Surface  Waves 

Measured  values  of  the  energy  as  a function  of  posi- 
tion of  the  receiver  transducer  for  ultrasonic  waves  re- 
flected from  a water-  aluminum  interface  at  the  angle 
Op  = 30.5”  are  shown  in  Fig.  9a.  The  solid  line  is  cal- 
culated from  Eq.  (8)  with  values  of  the  various  p.-iram- 
eters  as  follows:  ■h'  =-14.8  mm,  4>”  =-5  mm*,  a = 8 
mm,  A = 1,  B = O.OIG. 


i,min 


ng.  1.  CompirtJon  bflween  •<H»pl*cement*  and  - ♦'  (fttrt  derivative  of 
phase  fhilt  upon  rencrllon)  at  i function  of  Incident  angle  for  water-al- 
uminum Interface. 


i,mtn 


face. 


Fig.  9.  Energy  redistribution  In  a Gaussian  ultrasonic  beam  reflected  from: 
a)  a water— aluminum  Interface;  b)  a water— brass  interface;  = bp  in  both 
cases. 


Energy  redistribution  in  an  ultrasonic  Gaussian  beam 
from  a water— brass  Interface  is  shown  in  Fig.  9b.  The 
angle  of  Incidence  was  the  optimum  angle  for  surface 
wave  generation  Op  = 47.0°.  The  theoretical  curve  was 
calculated  from  Eq.  (8)  with  the  following  experimentally 
determined  parameters:  ♦'  = — 16.5  inm,  4>'  =-7  mm*, 
a = 8 mm,  A = 1,  and  B = O.OIG. 

Both  the  experimental  and  the  theoretical  curves  in 
Figs. 9a and 9b  indicate  that  an  incident  Gaussian  beam, 
upon  reflection  at  the  angle  of  optimum  generation  of 
surface  waves,  is  redistributed  into  multiple  beams  (two 
of  which  were  detcct.ahlc  over  the  distance  of  travcr.se  of 
the  receiver  transducer).  The  schllcren  photograph 
Figure  5b  indicates  that  there  are  actually  more  than  two 
reflected  beams. 

5.  SUMMARY 

It  has  been  suggested  in  this  paper  that  an  ultra.sonic 
beam  reflected  from  a liquid-  solid  interface  at  the  angle 
at  which  surface  waves  arc  generated  has  a different  en- 
ergy distribution  from  the  Incident  beam.  Therefore,  the 
concept  of  the  "beam  displacement"  expresses  only  a 
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gross  effect.  A(q>lyiDg  the  mathematical  analysis  which 
was  introduced  by  Brekhovskikh,  an  analytic  expression 
was  derived  for  the  reflected  amplitude  for  the  case  of 
an  Incident  Gaussian  beam.  The  theory  was  experimentally 
tested  for  water— aluminum  and  for  water— brass  inter- 
faces. The  transmitter  and  the  receiver  transducers  were 
speciaUy  designed  apodized  transducers  with  a Gaussian 
output.  Both  the  theoretical  and  the  experimental  results 
indicate  that  the  incident  Gaussian  beam  is  redistributed 
by  the  reflection  process  into  multiple  beams.  The  posi- 
tion and  the  relative  heights  of  these  beams  are  functions 
of  the  transducer  characteristics . They  depend  also  on 
the  liquid— solid  interface  through  the  first  and  second 
derivative  of  the  phase  shift  upon  reflection.  There  is 
fair  agreement  between  theory  and  experiment  concerning 
the  position  and  height  of  the  secondary  beam,  while  good 
agreement  is  obtained  for  the  position  of  the  principal 
beam. 

To  understand  the  "displacement"  phenomenon,  one 
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can  consider  the  finite  extension  of  the  beam  as  a super- 
position of  partial  plane  waves  whose  mutual  phase  rela- 
tions are  sensitive  functions  of  their  angle  of  reflection. 
The  individual  plane  wave  components  are  therefore  super- 
imposed with  mutual  phase  relations  which  differ  from 
those  of  the  incident  beam.  The  change  of  the  mutual  phase 
relations  of  the  plane  wave  components  is  interpreted  as 
the  physical  origin  of  the  "beam  displacement,"  or,  in 
fact,  the  energy  redistribution. 
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A Gaussian  distribution  of  amplitudes  in  an  ultrasonic  beam  reHected  from  liquid-solid  interfaces  is 
used  in  the  reexamination  of  the  concept  of  “beam  displacement”  which  occurs  at  the  angle  of 
excitation  of  surface  waves  on  the  interface.  Experimental  results  show  that  the  energy  of  the 
reflected  beam  is  redistributed  into  two  or  more  beams  at  (or  near)  this  angle.  The  theory  of 
Brekhovskikh  has  been  extended  to  include  both  a Gaussian  input  beam  and  the  second  derivative  of 
the  phase  shift  upon  reflection.  Reasonable  agreement  is  obtained  between  theory  and  experiment  for 
water-aluminum  and  water-brass  interfaces.  For  water-beryllium  and  water-stainless  steel  the 
agreement  is  fair. 

Subject  Classification:  35.24,  35.26;  20.30;  45.10. 


INTRODUCTION 

Consider  an  ultrasonic  beam  incident  on  a liquid- 
solid  interface.  For  a wide  range  of  angles  of  incidence 
this  problem  can  be  mathematically  described  by  solving 
the  boundary  value  problem  under  the  assumption  of 
energy  conservation.  At  the  critical  angles,  however, 
experimental  verification  of  the  calculated  results  is 
complicated  by  the  fact  that  the  reflected  beam  appears 
to  be  “displaced.  ” This  effect  is  especially  pronounced 
at  the  angle  for  excitation  of  surface  waves  along  the 
interface,  as  has  been  demonstrated  by  Schoch.  ‘ The 
"displacement  ” of  the  reflected  beam  was  attributed  to 
the  fact  that  a beam  does  not  behave  as  an  Infinitely  ex- 
tended plane  wave. 

The  situation  is  shown  in  Fig.  1.  An  ultrasonic  beam 
is  incident  upon  an  Interface  at  an  angle  6^  (approximate- 
ly the  angle  at  which  surface  waves  are  excited  on  the 
interface).  The  beam,  assumed  to  be  infinitely  extended 
in  a direction  normal  to  the  figure,  is  displaced  a dis- 
tance A upon  reflection.  In  developing  a theory  to  de- 
scribe this  effect,  one  defines  the  reflection  coefficient 
as  a complex  quantity  fle**'*’  and  includes  the  mutual 
phase  relation  of  partial  waves  in  the  reflected  beam  by 
expanding  the  expression  for  the  phase  shift  upon  re- 
flection 0(/))  into  a power  series  of  the  form* 

<<>(/>)  =<f> + <f>’(/>  - ®) + 2 -<*)*  + ••  • ( (1) 

where  <f)=<)(a),  (t>'={S<i>/Bp),.^,  = and 

P=k3ln6,  a=ksln6g,  k = 2‘it/\,  X = the  ultrasonic  wave- 
length in  the  liquid,  6g  is  the  angle  of  incidence,  and  6 
is  the  angle  of  incidence  of  a wavelet  in  the  incident 
beam.  By  carrying  out  this  analysis  in  the  first  ap- 
proximation, one  finds  that  the  displacement  is* 

A = -(8(j)/8/))^„  = -4,'  . (2) 

In  this  approximation,  the  beam  profile  of  the  re- 
flected beam  (l.e.,  the  energy  distribution  of  the  re- 
flected beam  measured  parallel  to  the  inierface)  is  the 
same  as  the  beam  profile  of  the  incident  beam.  Such 
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profiles  are  sketched  on  Fig.  I as  Fig.  1(a)  and  Fig. 

Kb). 

In  making  a comparison  between  theory  and  experi- 
ment, it  has  been  found*"*  that  there  is  nominal  agree- 
ment between  the  measured  “displacement”  and  that 
calculated  from  Eq.  1;  however,  the  experiment  is  com- 
plicated by  the  presence  of  unexpected  reflected  beams. 
For  this  reason,  it  is  necessary  to  carry  out  the  anal- 
ysis in  the  second  approximation  to  include  the  second- 
order  terms  in  Eq.  1.  In  this  higher  approximation,  the 
“beam  displacement  ” loses  its  meaning.  The  profile 
of  the  reflected  beam  is  not  the  same  as  that  of  the  in- 
cident beam;  i.e.,  there  is  an  “energy  redistribution” 
in  the  reflected  beam.  [Figure  1(b)  would  no  longer  be 
the  same  as  Fig.  1(a).  ] 

Brekhovskikh*  has  developed  a detailed  theory  for  the 
reflection  of  a bounded  ultrasonic  beam  at  a liquid- 
solid  interface.  From  this  theory  he  obtained  an  ex- 
pression for  the  amplitude  of  the  reflected  beam.  Bre- 
khovskikh also  pointed  out  that  the  structure  of  the  re- 
flected beam  is  in  general  different  from  that  of  the  in- 
cident beam,  i.e.,  that  there  is  an  “energy  redistri- 
bution. ” He  obtained  an  expression  for  the  amplitude 
distribution  in  the  reflected  beam  as  a function  of  the 
distance  along  the  interface,  and  solved  a special  case 
under  the  assumption  that  the  incident  beam  is  accurate- 
ly described  by  a rectangular  function. 

In  a comparison  between  theory  and  experiment,  we 
have  encountered  the  fact  that  a rectangular  function  dis- 
tribution in  the  incident  ultrasonic  beam  leads  to  dif- 
fraction effects  which  confuse  the  comparison.  For  this 
reason,  we  have  developed  a transducer  which  produces 
an  amplitude  distribution  in  the  Incident  beam  which  is 
accurately  described  by  a Gaussian  function.  * The 
Gaussian  function  amplitude  distribution  produces  a sin- 
gle well-defined  beam  without  diffraction  lobes.  The 
purpose  of  the  present  work  is  to  show  how  closely  the 
experimental  results  obtained  with  a Gaussian  distri' 
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Breazeale,  Adler,  and  Flax:  Rallection  trom  a liquid-solid  interface 


FIG.  1.  Reflection  of  an  ultrasonic  wave  at  a liquid-solid  in- 
terface. (a)  Incident  beam  cross  section,  (b)  First-order  rep- 
resentation of  reflected  beam  cross  section. 


button  of  amplitudes  in  the  incident  beam  agree  with 
those  calculated  from  the  appropriate  extension  of  the 
Brekhovskikh  theory. 
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I.  REFLECTION  OF  A GAUSSIAN  BEAM-THEORETICAL 
CALCULATION  OF  ENERGY  REDISTRIBUTION 

In  order  to  represent  the  beam  as  a superposition  of 
plane  waves,  Brekhovskikh  writes  the  field  in  the  plane 
z=l  (at  the  interface)  as  a Fourier  integral.  By  intro- 
ducing the  amplitude  reflection  coefficient  Klej),  which 
is  a function  of  the  angle  of  incidence  Sq  , the  general 
form  of  the  solution  to  the  wave  equation  which  de- 
scribes the  reflected  beam  is  written  as 

i/'r.fitr)  = ^ Y(eo)  exp{t(ax  + 1 V(^^-a^)]} 

X |exp^-t  ^u^^^F[x  + <P' + '^1^ u]du  , 

(3) 

where  I is  the  distance  between  the  transmitter  and  the 
liquid-solid  interface,  u = {x*<p' - £)/Vir^'  , £ is  the 


FIG.  2.  Experimental  arrange- 

(a)  ment.  (a)  Apparatus,  (b)  Dia- 

gram of  goniometer. 
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[11] 


A solution  of  Eq.  5 can  be  obtained  in  the  form  of 
error  functions  by  completing  the  square  and  writing  the 
two  exponential  functions  as  a single  exponential  func- 
tion whose  argument  is  squared.  The  result  is 


Interface 


(a)  Regular  Transducer 


where  we  have  used  the  definition  of  the  error  function 
2/ '/It  !o  exp(-  v^)  rft)  = 1.  The  energy  distribution  in  the 
reflected  beam  is 


This  is  basically  a Gaussian  fu.iction  whose  peak  has 
been  displaced  a distance  - 0'  from  the  ptrsition  x = 0, 
but  the  width  has  been  changed  by  the  retention  of  <&" 


|lnt»ffoce 


For  the  purposes  of  improving  agreemeiii  between 
theory  and  experiment,  one  can  also  define  limits  on 
the  integral*  as  follows : 


(bl  Apodized  Transducer 


FIG.  3,  Schlleren  photographs  of  transducer  outputs,  (a)  Rec- 
tangular function  amplitude  distribution,  (b)  Gaussian  amplitude 
distribution. 


where  a is  a measure  of  the  width  of  the  Gaussian  beam 
over  which  the  integration  is  taken.  In  this  case,  the 
reflected  beam  is  described  by 


variable  of  integration  in  the  original  Fourier  integral, 
and  F(x  + <l)'  + '/ )r<>"  u)  describes  the  incident  beam  func- 
tional dependence  on  the  amplitude  distribution  along  the 
interface.  (In  Eq.  3 the  higher  approximation.  Inclusion 
of  second-order  terms  in  Eq.  1,  has  been  made.) 


• AA«oiur«d 
~ Colculofcd  Goutiion 
Curve 


For  the  present  purposes.  Fix  + 0'  + u)  is  to  be 
a Gaussian  function.  Such  a function  does  not  change 
appreciably  over  a distance  of  one  wavelength  along  the 
propagation  direction.  This  satisfies  the  requirement 
Imposed  by  the  manner  in  which  Eq.  3 was  derived.’ 
The  Gaussian  function  describing  the  incident  beam  can 
be  written 


where  A and  B are  experimentally  determined  constants 
which  depend  on  the  characteristics  of  the  transducer. 
The  amplitude  of  the  reflected  beam  then  may  be  ob- 
tained by  substituting  Eq.  4 into  Eq.  3: 


FIG.  4.  Energy  distribution  in  the  ultrasonic  beam  emitted  by 
the  apodized  transducers. 
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i^'ref  1 


4^B4>"  + i/i 


exp  i a.v  + / -Jk^  - ■ 


2^<p  "(x  + 0 ')^  "I 

45^0  "^rj 


This  function  can  be  evaluated  by  computer  and  the  prod- 
uct 0*0  formed  to  evaluate  the  energy  distribution  as  a 
function  of  position  x.  Such  a function  is  no  longer  a 
simple  Gaussian  function. 

II.  REFLECTION  OF  A GAUSSIAN  BEAM- 
EXPERIMENTAL  MEASUREMENT  OF  ENERGY 
REDISTRIBUTION 

A.  Apparatus 

The  experimental  apparatus  is  shown  in  Fig.  2.  The 
electronic  system  is  a standard  pulse  arrangement, 
with  the  signal  from  the  receiver  transducer  being  dis- 
played on  an  oscilloscope.  Two  aspects  of  the  arrange- 
ment are  unique,  however:  the  goniometer  and  the 
transducers.  In  the  experimental  situation  both  the 
driver  and  the  receiver  transducer  are  mounted  on  the 
goniometer  and  this  unit  is  placed  in  a tank  of  water. 

/.  Goniometer  system 

A schematic  diagram  of  the  goniometer  is  shown  in 
Fig.  2(b).  The  system  is  designed  so  that  the  face  of  the 
specimen  (or  reflector)  forms  a vertical  plane  at  the 
center  of  a circle  formed  by  rotation  of  the  two  arms. 

A vernier  scale  on  each  arm  permits  the  angular  posi- 
tion of  each  arm  to  be  read  to  two  minutes  of  arc.  Both 
arms  have  transducer-holder  stages  mounted  on  riders 
that  allow  the  transducers  to  be  moved  up  to  50  cm  in 
a radial  direction.  In  addition,  the  transducer  on  the 
receiver  arm  can  be  moved  perpendicular  to  the  arm 
(or  radius).  This  way  accurate  scanning  of  the  re- 
flected beam  can  be  accomplished. 
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B(x  + 4>'}'f^'' 
Vb0"  + ?'/2  . 


,r  /„  .//  r?  B(x  + <P*) 'll  ,a\ 
-erf[iq>0  ^ 


I 

2.  Transducer 


Details  of  the  behavior  of  an  ultrasonic  wave  re- 
flected from  an  interface  depend  fundamentally  on  the 
energy  distribution  in  the  cross  section  of  the  incident 
ultrasonic  beam  (see  Eq.  3).  The  mathematical  inter- 
pretation of  the  results  is  facilitated  if  the  energy  dis- 
tribution in  the  Incident  beam  cross  section  can  be  de- 
scribed in  a simple  relation,  such  as  a Gaussian  func- 
tion. The  experimental  measurements  are  also  more 
straightforward  if  one  uses  a Gaussian  function  distri- 
bution in  the  incident  ultrasonic  beam.  This  can  be 
demonstrated  by  schlieren  photography. 

Figure  3(a)  is  a continuous -wave  schlieren  photo- 
graph of  a 4-MHz  ultrasonic  beam  from  a narrow  trans- 
ducer Incident  at  an  angle  of  35°  upon  a water-aluminum 
interface.  This  beam  has  a rectangular  function  distri- 
bution of  amplitudes  near  the  transducer.  The  interpre- 
tation problem  resulting  from  the  diffraction  lobes  is 
obvious.  Figure  3(b)  shows  an  ultrasonic  beam  of  ap- 
proximately the  same  width  and  frequency  with  a Gauss- 
ian amplitude  distribution  across  its  width.  This  beam, 
produced  by  a new  transducer,  ^ is  incident  on  the  same 
interface  at  the  same  angle  of  incidence.  In  addition  to 
the  absence  of  diffraction  lobes,  there  is  less  divergence 
in  the  Gaussian  ultrasonic  beam. 


The  Gaussian  distribution  was  obtained  by  using  a 
narrow  electrode  at  the  back  of  the  quartz  transducer 
and  grounding  the  entire  front  plane.  Fringing  of  the 
electrical  field  in  the  quartz  transducer  material  closely 
approximates  a Gaussian  function  if  the  ratio  of  the 
width  W of  the  back  electrode  to  the  thickness  T of  the 


FIG.  5.  "Displacement”  of  a Gaussian  beam  incident  upon  a 
water-aluminum  interface,  as  a function  of  the  angle  of  inci- 
dence. 


FIG.  6.  “Displacement”  of  a Gaussian  beam  incident  upon  a 
water-brass  interface  as  a function  of  the  angle  of  incidence. 
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FIG.  7.  “Displacement”  of  a Gaussian  beam  Incident  upon  a 
water-beryllium  interface  as  a function  of  the  angle  of  inci- 
dence. 


transducer  material  is  2<  W/T < 4. 

The  ultrasonic  beam  profile  emitted  by  these  trans- 
ducers could  be  obtained  by  using  the  goniometer  sys- 
tem with  the  sample  removed.  All  measurements  with 
the  goniometer  system  were  made  at  a frequency  of  2 
MHz.  A plot  of  the  energy  distribution  obtained  by  using 
the  new  transducers  both  as  transmitter  and  receiver 
is  given  in  Fig.  4.  (The  transmitter  is  approximately 
five  times  as  wide  as  the  receiver. ) The  experimental 
points  have  been  fit  by  a Gaussian  function  of  the  form 
where  A = 1 and  B = 0.016  cm'^.  The  agreement 
between  the  calculated  curve  and  the  experimental  points 
indicates  that  the  ultrasonic  beam  produced  by  the  new 
transducers  is  accurately  described  by  a Gaussian  func- 
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FIG.  8.  "Displacement”  of  a Gaussian  beam  incident  upon  a 
water-stainless  steel  interface  as  a function  of  the  angle  of  in- 
cidence. 


(c) 

FIG.  9.  Schlieren  photograph  of  Gaussian  ultrasonic  beam  re- 
flected from  a water-brass  interface,  (a)  95<e,j.  (b) 
energy  redistribution  is  evident,  (c) 


tion.  Since  secondary  diffraction  maxima  are  no  longer 
present,  one  could  accurately  describe  these  devices  as 
“apodized  tranducers. ” 

III.  RESULTS 

For  various  angles  of  incidence,  one  can  follow  the 
position  of  the  maximum  in  the  reflected  beam.  Within 
the  ^proximations  of  the  Schoch  theory,  a plot  of  the  pro- 
jection of  the  distance  moved  by  the  receiving  trans- 
ducer onto  the  plane  of  the  face  of  the  reflection  would 
be  the  “displacement.”  Plots  of  this  quantity  for  alumi- 
num, brass,  beryllium,  and  stainless  steel  are  given  in 
Figs.  5-8.  Superimposed  on  the  experimental  data  are 
plots  of  -<!>',  as  indicated*  in  Eq.  2.  The  critical  angle 
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for  the  longitudinal  wave  in  the  solid  for  the 

shear  wave  in  the  solid  (8„),  as  well  as  the  optimum  angle 
for  excitation  of  surface  waves  on  the  interface  are 
indicated.  It  is  characteristic  that  the  values  of  - 0 ' at 
the  angle  for  excitation  of  surface  waves  6g  are  consis- 
tently greater  than  the  measured  "displacement.  ” This 
results  from  the  fact  that,  at  6^,  there  is  in  fact  an 
energy  redistribution,  rather  than  a simple  “displace- 
ment” of  the  reflected  beam. 

Careful  schlieren  photography  has  shown  us  that  even 
with  the  apodized  transducers,  which  produce  a Gauss- 
ian distribution  of  energy  in  the  incident  beam,  the  beam 
is  reflected  as  multiple  beams  when  the  angle  of  inci- 
dence is  that  at  which  surface  waves  are  excited  on  the 
interface.  This  is  shown  in  Fig.  9(b)  for  a water- 
brass  interface,  in  which  the  null  strip  described  by 
Neubauer’  is  visible.  Figure  9(a)  and  9(c)  are  schlieren 
photographs  for  angles  of  incidence  respectively  smaller 
and  larger  than  this  angle. 

Scans  of  the  reflected  beams  with  the  goniometer  sys- 
tem for  angles  respectively  smaller  and  larger  than  the 
angle  for  surface  wave  excitation  gave  undisplaced 
Gaussian  distributions.  The  results  for  a brass  reflec- 
tor shown  in  Fig.  10  are  typical.  The  curves  are  Gauss- 


nSTANCC ALONG  NTERFACE  (rmi) 


FIG.  11.  Energy  redistribution  in  an  ultrasonic  beam  reflected 
from  a water-aluminum  Interface  at  5*.  Distance 

is  measured  from  the  center  of  the  Incident  beam. 


ian  functions  of  the  same  form  as  that  for  the  incident 
beam  shown  in  Fig.  4.  An  incident  angle  of  10°  is  well 
below  the  longitudinal  critical  angle  for  brass;  an  angle 
at  54°  is  well  above  the  angle  for  excitation  of  surface 
waves.  A straightforward  calculation  from  the  simpli- 
fied theory  gives  an  energy  reflection  coefficient  of  80% 
at  10°  and  100%  at  54°.  These  values  are  consistent 
with  the  maximum  values  of  the  curves  in  Fig.  10.  The 
other  reflectors  produced  similar  results. 

Scans  of  the  reflected  beam  with  the  incident  beam 
striking  the  surface  at  the  angle  6^  allowed  a measure  of 
the  energy  redistribution.  These  results  can  be  com- 
pared directly  with  the  theoretical  predictions  from 
either  Eq.  7 or  8,  if  one  separates  out  the  effect  of  the 
"specular”  beam  described  by  Neubauer.’ 

Figures  11-14  give  the  experimental  data  points  obtained 
by  squaring  the  magnitudes  measured  across  the  re- 
flected beam,  correcting  for  the  fact  that  the  receiver 
moves  perpendicular  to  the  radius  rather  than  parallel 
to  the  Interface.  Special  attention  was  paid  to  making 
the  receiver  angle  exactly  equal  to  the  angle  of  inci- 
dence. Because  of  this,  the  maximum  is  consistently 
"displaced”  a smaller  distance  than  expected  from  the 
value  of  0 ' at  the  angle  , or,  indeed,  from  the  ex- 
perimental data  in  Figs.  5-8. 
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FIG.  13.  Energy  redistribution  In  an  ultrasonic  beam  reflected 
from  a water-beryllium  interface  at  8g  = Aj,°12*. 
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Superimposed  on  the  experimental  data  in  Figs.  11-14 
are  theoretical  curves  calculated  from  Eq.  8.  For  ref- 
erence, values  of  the  “displacement  ” Aj  predicted  by 
the  Schoch  theory  are  indicated.  The  fact  that  the  peaks 
in  the  experimental  data  are  consistently  wider  than  the 
theoretical  curves  is  probably  associated  with  the  fact 
that  the  receiving  transducer  does  not  have  the  ideal 
infinitesimal  width.  Values  of  the  parameters  <^)',  </>”, 
a,  and  B used  in  the  calculation  for  each  material  are 
listed  in  Table  I.  The  values  of  <p ' were  determined  by 
the  positions  of  the  maxima  in  the  experimental  data. 

The  value  of  0"  is  not  readily  measured;  thus,  a value 
was  selected  which  produced  best  agreement  between 
theory  and  experiment.  The  value  of  B was  determined 
by  projecting  the  Gaussian  function  shown  in  Fig.  4 
onto  the  interface.  The  value  of  a came  from  the  follow- 
ing considerations. 

The  photograph  of  Fig.  9(b),  as  well  as  the  experi- 
mental data  of  Figs.  11-14,  show  secondary  beams  in 
addition  to  the  main  beams.  Such  secondary  beams  are 
not  predicted  by  the  theory— even  including  the  term  in 
0”— unless  the  Gaussian  function  is  truncated.  A finite 
value  of  a is  a measure  of  the  position  at  which  this 
truncation  occurs.  The  value  a = 10  gave  the  secondary 
peak  at  34  mm  in  Fig.  11  for  the  aluminum  data  and 
improved  agreement  near  - 5 mm.  This  value  was 
maintained  for  the  other  curves  to  determine  the  extent 
to  which  the  secondary  peaks  would  be  predicted  by  theory 
for  these  materials  as  well. 


IV.  SUMMARY  AND  CONCLUSIONS 

The  behavior  of  an  ultrasonic  wave  at  a liquid-solid 
interface  is  reasonably  well  described  by  a solution  of 
the  straightforward  boundary  value  problem  at  most 
angles  of  incidence.  At  the  angle  of  incidence  lor  which 
the  surface  waves  are  excited,  this  description  is  in- 
adequate. Schoch’s  introduction  of  the  idea  of  a “dis- 
placement, ” while  an  improvement,  also  proves  to  be 
inadequate.  At  this  angle  one  observes  that  the  reflected 


FIG.  14.  Energy  redletributlon  In  an  ultrasonic  beam  reflected 
from  a water-stainless  steel  Interface  at  8g  = ° 39. 8*. 


TABLE  I.  Parameters  used  in  calculating  the 
theoretical  curves. 


a 

B 

Aluminum 

-9 

-33 

10 

0.012 

Brass 

-12.4 

-43 

10 

0.0076 

Beryllium 

-2 

-13 

10 

0.016 

Stainless  steel 

-7 

-18 

10 

0.012 

beam  has  a different  energy  distribution  from  the  inci- 
dent beam.®  Therefore,  the  concept  of  “beam  displace- 
ment ” expresses  only  a gross  effect. 

The  existence  of  a transducer  which  produces  a truly 
Gaussian  distribution  of  amplitude  has  made  it  possible 
to  make  a direct  comparison  between  experimental  re- 
sults from  four  different  interfaces  and  a theory  based 
on  the  approach  of  Brekhovskikh.  Both  the  experimental 
and  the  theoretical  results  indicate  that  the  incident 
Gaussian  beam  is  redistributed  in  the  reflection  process 
into  multiple  beams.  The  positions  and  relative  in- 
tensities of  these  beams  are  functions  of  the  transducer 
characteristics.  They  depend  on  the  characteristics 
of  the  liquid-solid  interface  through  the  first  and  second 
derivative  of  the  phase  shift  upon  reflection.  The  sec- 
ondary beams  also  seem  to  be  associated  with  truncation 
of  the  incident  beam.  There  is  reasonable  agreement 
between  theory  and  experiment  for  the  aluminum  and  the 
brass  interfaces,  while  the  agreement  for  the  beryllium 
and  the  stainless  steel  interfaces  is  less  good.  In  all 
cases,  the  experimental  value  of  the  “displacement  ” is 
different  from  that  predicted  by  Schoch.'  This  dif- 
ference is  quite  great-a  factor  of  26-for  beryllium. 
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The  behavior  of  a Gaussian  ultrasonic  beam  incident  on  a liquid-solid  interface  at  the  Rayleigh  angle,  the 
angle  at  which  surface  waves  are  excited  on  the  interface,  has  been  studied  in  some  detail.  The  reflected 
beam  is  displaced  in  the  manner  predicted  by  Schoch;  however,  the  “Schoch  displacement”  in  general  is 
too  large.  Good  agreement  is  obtained  between  experimental  results  and  the  theory  of  Bertoni  and  Tamir, 
which  assumes  that  the  incident  beam  couples  resonantly  into  a leaky  surface  wave  at  the  Rayleigh  angle 
and  that  the  energy  reradiated  from  this  leaky  surface  wave  interferes  with  specularly  reflected  energy.  The 
propagation  distance  of  the  ultrasonic  beam  is  explicitly  included  in  describing  the  ultrasonic  wave 
reflection  at  the  Rayleigh  angle. 

PACS  numbers:  68.25. -|-j,  43.35. -(-d,  68.30. -t-z 


I.  INTRODUCTION 

A wave  incident  on  a plane  interface  can  experience 
total  reflection  under  certain  conditions . At  the  smallest 
angle  at  which  total  reflection  occurs,  one  part  of  the 
reflected  (or  refracted)  beam  propagates  tangentially 
along  the  surface.  Interaction  between  this  tangential 
component  and  other  (reflected  or  refracted)  components 
leads  to  a lateral  displacement  of  the  reflected  beam. 
This  interaction  takes  place  because  the  wave  has 
penetrated  into  the  second  medium — even  under  con- 
ditions of  total  reflection. 

Even  though  these  observations  are  consistent  with 
electromagnetic  theory,  and  even  though  Newton*  in 
1704  suspected  that  a light  wave,  upon  total  internal 
reflection  at  a glass-air  interface,  penetrates  into  the 
air,  a direct  demonstration  of  the  effect  did  not  appear 
until  Goos  and  Hanchen*  modified  Newton’s  original 
experiment  by  bringing  a silver  layer  into  contact  with 
a portion  of  the  glass  surface  at  which  total  reflection 
took  place.  Lateral  displacement  of  the  totally  reflected 
beam  was  demonstrated  by  direct  comparison  with  a 
portion  of  the  same  beam  reflected  from  the  silvered 
area  in  which  negligible  penetration  occurred.  The  effect 
is  greatest  at  the  critical  angle,  and  becomes  less 
pronounced  for  larger  angles  of  total  reflection.  An 
extensive  survey  of  this  subject  is  given  by  Lotsch,® 
who  also  mentions  the  acoustical  analogue. 

An  ultrasonic  wave  incident  upon  a liquid-solid  inter- 
face exhibits  two  critical  angles — one  for  the  longitudinal 
wave  in  the  solid  and  one  for  the  shear  waves . For 
angles  much  larger  than  the  critical  angle  for  shear 
waves , total  reflection  also  is  observed . Near  the 
shear  wave  critical  angle,  the  reflected  beam  is  dis- 
placed in  a manner  analogous  to  the  optical  beam.  But 
it  turns  out  that  a different  mechanism  is  responsible 
for  most  observed  displacements  of  ultrasonic  waves. 

By  using  schlieren  photography,  Schoch*  wss  able  to 
demonstrate  a displacement  of  the  reflected  beam  by 
photographing  the  ultrasonic  beams  in  the  liquid.  As 
shown  in  Fig.  1,  the  reflected  beam  is  displaced  a 
distance  A laterally  along  the  interface;  however,  the 
angle  of  maximum  displacement  is  somewhat  greater 
than  the  critical  angle  for  shear  waves.  Schoch  also 


developed  a theory  which  predicted  the  lateral  dis- 
placement of  the  ultrasonic  wave  after  reflection.  The 
“Schoch  displacement”  is  attributed  to  the  fact  that  the 
physically  realizable  ultrasonic  beam  does  not  behave 
as  an  infinitely  extended  plane  wave  as  is  usually  as- 
sumed in  solution  of  boundary-value  problems.  Schoch 
included  the  mutual  phase  relationship  of  partial  waves 
in  the  reflected  beam  by  expanding  the  expression  for  the 
phase  shift  upon  reflection  into  a power  series . The 
first  derivative  of  the  phase  shift  is  identified  as  the 
“Schoch  displacement” . 

Brekhovskikh*  extended  Schoch’s  theory,  and  by  in- 
cluding the  second  derivative  of  the  phase  shift  he  de- 
rived an  expression  for  the  amplitude  distribution  of 
the  reflected  beam.  This  theory  subsequently  was  modi- 
fied to  account  for  a Gaussian  amplitude  distribution  in 
the  incident  biam,  and  an  apodized  transducer  which 
produced  a Gaussian  amplitude  distribution*  in  the  inci- 
dent beam  was  used  in  a comparison  between  theory 


FIG.  1.  Schlieren  photograph  of  an  ultrasonic  beam  incident  on 
a liquid-solid  (aluminund  interface  at  the  Rayleigh  angle.  The 
reflected  beam,  as  indicated,  is  split  into  two  components:  a 
specular  beam  and  a beam  displaced  a distance  laterally  down 
the  Interface.  Secondary  beams  are  visible  at  greater  distances. 
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and  experiment.’  A closed-form  analytical  solution 
improved  the  theory  and  simplified  the  calculation,® 
Nominal  agreement  between  theory  and  experiment  was 
obtained  for  water-aluminum  and  water-brass  inter- 
faces, and  poor  agreement  was  obtained  for  water- 
beryllium  and  water— stainless -steel  interfaces.  In  spite 
of  nominal  agreement  between  theoi-y  and  experiment, 
these  theories  were  not  totally  satisfactory  since  the 
physical  mechanism  responsible  for  the  displacement 
did  not  enter  explicitly.  Furthermore,  they  accounted 
for  the  displaced  beans  only.  A number  of  theorists®'®"” 
and  experimenters’®"*®  have  correctly  associated  the 
displaced  beam  with  the  excitation  of  surface  waves  on 
the  interface,  have  pointed  out  that  Lord  Rayleigh 
studied  surface  waves  on  a free  solid  interface,  and 
have  labelled  the  angle  at  which  the  displacement  is  a 
maximum  9^  to  distinguish  it  from  the  critical  angle 
for  shear  waves  in  the  solid. 

Experimental  studies  of  the  reflection  of  ultrasonic 
waves  from  such  solids  as  beryllium  and  AljOg  revealed 
a second  reflected  beam  which  is  not  displaced,  but  is 
specularly  reflected.  This  beam  also  is  visible  in  Fig, 

1 , It  Is  separated  from  the  displaced  beam  by  a “null 
strip’’  since  the  two  beams  are  out  of  phase,*® 

Bertoni  and  Tamir*®  introduced  a new  theory  to  ex- 
plain the  acoustical  Rayleigh-angle  problem.  This 
theory  is  based  on  earlier  work*®  to  describe  the  Goos- 
Hanchen  effect  in  optics.  Their  postulated  mechanisms 
include  a specular  wave  and  a reradiating  surface  wave, 
which  they  call  a “leaky  wave”  because  it  leaks  energy 
back  into  the  liquid  as  the  displaced  beam.  They  present 
a solution  for  the  reflected  field  which  is  strictly  valid 


only  if  both  source  and  point  of  observation  are  at  the 
interface.  In  the  present  work  we  point  out  that  it  is 
necessary  to  include  the  ultrasonic  wave  propagation 
distance  (from  the  emitting  transducer  to  the  interface 
to  the  receiving  transducer)  in  order  to  account  for  the 
effect  of  diffraction.  Although  it  is  not  necessary  to 
include  the  propagation  distance  in  the  analogous  optical 
problem,  Horowitz  and  Tamir”  recognized  that  it 
could  be  included.  Since  the  dependence  of  this  distance 
is  omitted  from  the  theory  of  Bertoni  and  Tamir,*®  we 
show  how  it  enters.  A detailed  comparison  is  made  of 
the  theoretically  predicted  energy  distribution  with  ex- 
perimental results  obtained  with  an  ultrasonic  gonio- 
meter using  apodized  transducers, 

II.  THEORY 

An  incident  bounded  beam  is  represented*®"*’  by  a 
Fourier  transform  pair  as 

= Y(^,)exp[f(/?^Ar+fc^)]rf/f^  (1) 

and 

Vi^^ix,0)exp(,-ik^)dx,  (2) 

where  and  are  the  x and  z components  of  the  wave 
number  k=  2t/\.  ri„j(x,0)  is  the  amplitude  distribution 
of  the  incident  beam  at  the  liquid-solid  interface.  A 
solution  for  the  reflected  wave  can  also  be  written  in 
the  form  of  a Fourier  integral  as 

exp[«(fejX  - k^z)]dk^. 

(3) 


Here  R ik,)  is  the  reflection  coefficient  given  as 


where  and  are  the  wave  numbers  corresponding  to 
the  longitudinal  (compressional)  and  shear  (transverse) 
waves,  respectively,  in  the  solid  and  p is  the  density 
ratio  of  solid  to  liquid.  R(fe,)  displays  a complex  char- 
acter which  implies  that  k,  can  be  generalized  to  a 
complex  variable  and  R(kJ  can  be  represented  in  the 
complex  plane.  Bertoni  and  Tamir  recognized  the  com- 
plex character  of  and  investigated  its  singular  points . 
The  integral  in  Eq.  (3)  is  evaluated  along  the  real  axis, 
so  singular  points  on  and  near  the  axis,  particularly 
those  whose  real  parts  lie  between  and  k,  will 
strongly  affect  its  value.  The  pole  near  the  real  axis  in 
this  interval  is  located  at 

fe^=fesine^ +i«,  (5) 

where  6^  is  near  the  Rayleigh  angle.  This  pole  repre- 
sents a resonance  of  the  system  in  the  form  of  a prop- 
agating wave:  exp{ik^x)  = exp(ik  sinO^x  - ax) . This  wave 
is  similar  to  Rayleigh  waves,  except  that  the  presence 
of  the  imaginary  part  a in  Eq.  (5)  causes  the  wave  to 


attenuate  as  it  progresses  along  the  x direction.  Energy 
is  continuously  “leaking”  away  from  the  solid  region  to 
the  liquid  region.  This  leaky  wave,  then,  is  a modified 
form  of  the  Rayleigh  wave. 

From  the  experimental  results  we  wish  to  report, 
one  cannot  reliably  evaluate  a because  of  diffraction 
effects.  We  therefore  follow  Bertoni  and  Tamir’®  and 
introduce  the  relationship  between  a and  the  Schoch 
displacement 

2o . . 

• ~ (fesinS,  - /tsinfl^)®  + a®  ’ ' ' 

where  S,  is  the  angle  of  incidence.  As  will  be  shown, 

Aj  can  be  calculated  from  measurable  physical  quanti- 
ties, If  0,  = 6^,  then 

or  = 2/ A,.  (7) 

Bertoni  and  Tamir*®  introduced  a simpler  function  for 
fl(fr,) — in  order  to  carry  out  an  analytical  integration 
for  Eq.  (3) — in  the  form 
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FIG.  2.  Coordinate  systems.  Typical  intensity  distributions 
are  superimposed  for  incident  imd  reflected  beams. 


where  k*  is  the  complex  conjugate  of  k^.  Equation  (8) 
results  from  a Laurent  expansion  around  the  pole 
and  therefore  is  a good  approximation  only  in  the  region 
near  the  pole,  the  region  of  interest  in  our  experiment. 

then  is  written  in  the  form 


R(fe,)  = Ro(i^,)  + R, ()?,). 

(9) 

where 

f<i  Kf, 

(10) 

k^~k,  k,-k, 

(11) 

where  k^  is  a specific  value  of  corresponding  to  a 
particular  angle  of  incidence,  and  kg  is  identical  to  k* 
if  the  tosses  are  negligible.  In  our  comparison  with 


experiment  we  assume  losses  are  negligible  cor- 
responds to  the  geometric  reflection  and  R,  accounts 
for  the  leaky-wave  contribution  The  reflected  field  can 
be  calculated  if  the  analytic  form  of  the  incident  beam 
is  given.  For  this  experiment  a specially  designed® 
Gaussian  transducer  was  used,  and  the  incident  beam 
is  described  as  a Gaussian  function: 


exp[-  (y/ Wq)^  + ik{x\ 

■ir‘‘'’u>ocos0j 


(12) 


where  is  the  half -width  of  the  area  on  the  interface 
illuminated  by  the  incident  beam.  Substituting  Eq.  (12) 
into  Eq . (2)  as  an  initial  condition  yields 


Vik  ) = exp[-(fe, 

" cose, 

which,  along  with  Eq.  (9),  is  inserted  into  Eq.  (3)  to 
find  the  reflected  field.  The  resulting  reflected  field 
is  written  in  two  terms  as 


= + (14) 

where 

V (x  ^ (ln-ko\  r-  exp[-(fe:,-fe,)"0„)^l 

FoU.zi  27rU, -frJJ-  cose, 


and 


X exp[f  (fe^^r  - k^z)  dfe. 


(15) 


V^ljc,z)  = 


1 (h-ko\  f”  /kr-  fe,\  exp[-  (k,  - ki)^(iwof] 

cose, 


X exp[  i {k^x  -k^)]dk^. 


(16) 


The  field  represents  a specular  or  geometrical 
acoustic  reflection,  the  type  which  occurs  with  varying 
amplitude  at  all  angles . Its  integrand  is  unaffected  by 
the  poles  at  The  field  F,  represents  the  leaky-wave 
contribution,  whose  integrand  is  strongly  affected  by 
the  pole  at 

Bertoni  and  Tamir*’  evaluated  the  integrals  in  Eqs. 
(15)  and  (16)  for  ^ =0,  i.e.,  at  the  interface.  For  2^0, 
which  describes  our  experimental  situation,  the  inte- 
grals can  be  evaluated  only  if  the  relationship  between 
and  can  be  found.  For  the  electromagnetic  case, 
Tamir  and  Bertoni,^®  following  the  approach  of 
Horowitz  and  Tamir,”  approximated  k,  by  the  first 
three  terms  of  a Fresnel  expansion 

fe  =fecose-(fe  -fesine)tane-(^*~-^^r^^  ^ . (17) 

* \ 2*co8®e  / 

(The  physical  assumption  in  arriving  at  this  expression 
is  that  only  those  plane  wavelets  are  considered  whose 
propagation  direction  is  essentially  the  same  as  the 
direction  of  the  central  beam.) 


TABLE  I.  Properties  of  solids  used  In  calculation  of  Schoch  displacement. 


Material 


P 


(xio*  cm/scc) 


I's 

(xio®  cm/ sec) 


Aluminum  2.71  6.40  3.106 

Stainless  steel  7.85  5.80  3.130 

Aluminum  oxide  4.0  10.70  6.  .360 


(xio^  cm/ sec) 

2.848 

2.827 

5.690 


imm) 


15.87 

43.10 

88.13 
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refers  to  a lateral  shift  of  the  geometrical  center  of 
the  reflected  beam  and  was  predicted  by  Schoch  as 


FIG.  3.  Goniometer  used  in  making  intensity  measurements. 


By  substituting  Eq.  (17)  into  Eqs,  (15)  and  (16)  one 
can  solve  the  integral  exactly  (see  Appendix  of  Ref.  17). 
The  relationships  between  the  fixed  and  reflected  co- 
ordinate systems  shown  in  Fig.  2 are 

x^=xcos0-(h-z)sin6,  (18) 


z^=xsine+  {h-z)co88. 


These  relationships  and  the  results  of  Tamir  and 
Bertoni**  are  used  to  evaluate  the  expressions  for  Vg 
and  Fj  in  the  reflected  coordinate  system: 


exp[-  (y,/w,)^  + tfeix,! 
7r*^’ic,cos0, 


(19) 


V,U,«,)  = 2Fo(x^,z,)  ^1  - 2-^^exp(>'’)erfc(y)^,  (20) 
where 

^ +f  — (sine, -sine_)u»  (21) 

and 

1 - (z,-x,tane,)j  (22) 

Here  erfc  is  the  complementary  error  function,  w is 
the  half -width  of  the  incident  beam,  and  is  the 
Schoch  displacement  which  is  characteristic  of  the 
interface.  Since  these  expressions  are  to  be  compared 
with  experiments  performed  at  the  Rayleigh  angle 
(l.e. , for  Of  = 9g),  we  have  used  Eq.  (7)  to  relate  « and 

The  total  reflected  field  at  any  point  is  given  by  the 
(vector)  sum  of  two  complex  numbers  F,  and  F^.  The 
field  amplitude  is  equal  to  the  modulus  of  the  sum;  the 
relative  phase  is  computed  from  the  quotient  of  the  real 
and  imaginary  parts  of  Vg  + Vf,  Both  field  amplitude 
and  relative  phase  were  calculated  by  use  of  FORTRAN. 
Energy  is  proportional  to  the  square  of  the  amplitude. 

One  of  the  important  parameters  which  influences  the 
distribution  of  the  reflected  amplitude  is  the  so-called 
Schoch  displacement,  A,  [see  Eqs.  (20)  and  (21)].  A, 


2x  /r(y-s)\"^’‘ri  + 6sHl-q)-2s(3-2q)1 

■n  ^\s(s  - 1) ) L s - q J ’ 

(23) 


where 


(24) 


F is  the  velocity  in  the  liquid,  F^  and  F,  are  the  shear 
and  longitudinal  velocities  in  the  solid,  F„  is  the 
Rayleigh  velocity,  p is  the  ratio  of  the  densities  of  the 
solid  and  the  liquid,  respectively,  and  x is  the  wave- 
length of  the  sound  in  the  liquid.  From  physically  mea- 
surable quantities  the  Schoch  displacement  can  be 
calculated. 

Three  materials  were  used  in  this  experiment— alu- 
minum oxide,  aluminum,  and  stainless  steel.  All  were 
immersed  in  water.  The  calculated  Schoch  displace- 
ments at  2 MHz  for  these  materials  are  listed  in  Table 
I.  These  values  of  the  Schoch  displacement  were  used 
in  the  calculations  of  the  reflected  field  from  Eq.  (14). 
As  will  be  seen  in  the  data,  the  Schoch  displacement 
is  a convenient  mathematical  concept,  but  it  does  not 
accurately  represent  the  lateral  shift  of  the  reflected 
beam. 


FIG.  4.  Amplitude  and  phase  profile  of  Incident  beam. 
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III.  EXPERIMENTAL  APPARATUS 

A.  Goniometer 

The  experimental  apparatus  consists  of  a two-crystal 
system.  Pulsed  ultrasonic  waves  are  produced  by  a 
standard  gated  amplifier  and  are  displayed  on  an  oscil- 
loscope. The  transducers  are  mounted  on  a goniometer 
in  such  a manner  that  the  ultrasonic  beam,  reflected 
from  the  metal  interface  of  interest,  is  scanned  by  the 
movable  receiver  transducer.  Figure  3 shows  the 
goniometer  removed  from  its  water  tank.  For  operation 
it  is  immersed  in  water  up  to  a point  about  halfway  be- 
tween the  top  of  the  specimen  holder  and  the  base  of 
the  column  on  which  the  arms  mount . The  left  arm 
carries  the  transmitter  and  the  right  arm  carries  the 
receiver  (the  slide  moves  the  receiver  perpendicular  to 
the  arm).  The  transducer  carriages  can  be  moved 
radially  along  the  arms  and  clamped  at  radii  between 
5 and  50  cm. 

The  specimen  holder  accommodates  rectangular 
blocks  (2x4  in.)  up  to  1 in.  thick  and  right-circular 
cylinders  (3  in.  in  diameter)  up  to  a few  inches  long. 
Either  specimen  type  is  restrained  so  that  its  front 
(reflecting)  face  coincides  with  the  rotation  axis  of  the 
arms.  Figure  1 also  shows  a cylindrical  aluminum  oxide 
specimen  in  place.  The  aluminum  oxide  cylinder  is 
2.5  in.  in  diameter,  so  an  aluminum  ring  is  used  to 
adapt  it  to  the  holder. 

B.  Transducers 

For  the  present  investigation,  the  transducers  were 
l-in.-diam  x-cut  quartz  crystals  plated  with  conducting 
material  (gold  or  silver)  on  one  side.  The  plated  side 
was  used  at  the  front,  or  water,  side.  The  rear  elec- 
trode was  a metal  strip  aligned  to  be  symmetric  about 
a diameter  of  the  crystal,  parallel  to  the  goniometer 
axis . This  electrode  arrangement  produced  an  ultra- 
sonic beam  with  a Gaussian  profile.*  A 6.35-mm  elec- 
trode was  used  in  the  transmitter  to  produce  a narrow 
beam  with  good  signal  amplitudes  at  the  receiver.  The 
receiver  had  a 3.68-mm  electrode  which  consistently 
produced  narrow  beam  profiles.  Figure  4 is  the  ampli- 
tude and  phase  profile  made  with  these  transducers  at 
2 MHz.  The  amplitude  profile  is  fitted  with  a Gaussian 
curve  of  the  form  A exp[-  (x  - n)*/b^].  This  Gaussian 
amplitude  variation  across  the  wavefront  makes  possible 
a detailed  comparison  between  theory  and  experiment, 
since  it  is  the  case  that  was  solved  analytically. 

C.  Amplitude  maanirement 

An  Arenberg  unit  was  used  in  its  pulsed-amplifier 
mode.  It  functions  as  a gate,  turning  the  input  from  the 
external  variable  frequency  oscillator  on  and  off  and 
feeding  it  Into  the  amplifier  section.  The  output  is  a 
tone  burst,  l.e. , a burst  of  2 MHz  rf  about  16  psec 
long.  Amplitudes  were  measured  by  reading  the  voltage 
developed  across  the  receiver  crystals  from  an 
oscilloscope. 

D.  Phate  maaturamant 

Several  potential  applications  motivated  efforts  to 
measure  phase— (1)  goniometer  alignment,  (2)  tracing 
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wavefront  contours  in  transducer-beam  profiles,  and 
(3)  tracing  wavefront  contours  in  Rayleigh-angle  reflec- 
tions. A technique  devised  for  detecting  a constant  phase 
led  to  improvements  which  allowed  approximate,  but  re- 
peatable, measurements  of  relative  phase.  Moving  the 
oscilloscope  trigger  from  the  Arenberg  trigger  output 
to  the  transmitter  cable  causes  the  transmitter  pulse 
to  trigger  the  oscilloscope  sweep  and  allows  relative 
phase  measurements. 

The  Arenberg  trigger  output  is  coincident  with  its 
gate,  but  these  signals  are  both  timed  by  a rate  gen- 
erator in  the  Arenberg  which  is  not  synchronized  with 
the  external  rf  input  (from  the  VFO).  Therefore,  the 
gate  and  trigger  pulse  do  not  occur  at  the  same  point 
(or  phase)  on  the  rf  waveform  each  time  the  gate  opens . 
If  the  oscilloscope  sweep  is  triggered  by  the  Arenberg 
gate,  the  envelope  of  the  pulse  remains  stationary  on 
the  CRT  trace,  but  the  rf  waveform  jitters  back  and 
forth  inside  the  envelope.  When  the  oscilloscope  sweep 
is  triggered  directly  by  the  transmitter  pulse,  it  syn- 
chronizes with  the  rf  waveform,  and  the  waveform  re- 
mains stationary  on  the  CRT  while  its  envelope  jitters. 

In  this  latter  condition,  the  displayed  signal  from  the 
receiver  is  locked  in  phase  with  the  transmitter  driving 
signal.  All  of  the  pulses  are  coherent  with  one  another 
because  they  originate  from  the  same  cw  source,  and 
the  oscilloscope  trigger  circuits  fix  the  phase  point  on 
the  rf  waveform  at  which  the  sweep  starts.  The  only 
remaining  phase  changes  result  from  changes  in  the 
propagation  delay  of  the  acoustic  signal  along  its  water 
path,  and  these  are  observable  on  the  CRT  as  move- 
ment of  the  entire  waveform  across  the  screen. 

The  slide  motion  of  the  receiver  confines  it  to  a plane 
in  the  propagation  path  of  the  acoustic  wavefronts . The 
wavefronts  are  typically  curved,  so  there  will  be  a point 
of  tangency  with  the  receiver  plane  as  the  front  ad- 
vances . The  signal  detected  by  the  receiver  at  the  point 
of  tangency  occurs  earlier  in  time  than  signals  detected 
elsewhere  in  the  receiver  plane,  and  it  appears  farthest 
left  on  the  CRT  trace.  This  point  was  designated  the 
“maximum  phase  excursion",  and  was  defined  as  zero 


FIG,  5.  Reflected  intensity  distribution  for  a water— aluminum 
oxide  interface. 
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FIG.  6.  Reflected  intensity  distribution  for  a water— stainless 
steel  interface  at  a,  = 800  mm. 


relative  phase  for  all  experimental  measurements. 
Changes  in  relative  phase  were  recorded  in  90°  in- 
crements according  to  the  distance  the  waveform  moved 
to  the  right  on  the  CRT.  All  relative  phases  were 
recorded  as  negative  with  respect  to  the  maximum  phase 
excursion. 

IV.  EXPERIMENTAL  PROCEDURE 

Reports  of  other  experimenters  established  a range 
of  values  where  the  Rayleigh  angle  for  a material  was 
expected  to  lie.  Starting  at  angles  1—2°  below  the  ex- 
pected range,  the  transmitter  and  receiver  arms  were 
set  at  equal  angles.  Then  the  receiver  was  scanned 
through  the  full  range  of  the  cross  slide,  and  the  loca- 
tions and  values  of  all  local  maxima  and  minima  in  the 
amplitude  were  noted.  The  slide  was  set  to  a position 
corresponding  approximately  to  the  center  of  the  inci- 
dent beam.  The  oscilloscope  was  set  to  indicate  zero 
relative  phase  at  this  position,  and  relative  phase  was 
observed  as  the  slide  was  scanned  away  to  the  right  and 
left.  Locations  of  large  or  sudden  changes  in  phase 
were  noted. 

The  arms  were  repositioned  and  the  scanning  proce- 
dure repeated  at  intervals  of  15'  of  arc.  As  the 
Rayleigh  angle  was  approached,  the  maxima  and  minima 
occasionally  shifted  position  and  relative  value  slightly, 
and  one  minimum  would  become  an  absolute  minimum, 
i.e. , its  amplitude  would  decline  to  a value  less  than 
those  of  the  other  minima.  If  the  Rayleigh  angle  was 
passed  over,  the  amplitude  at  the  absolute  minimum 
would  pass  through  a minimum  of  its  own  and  start  to 
increase  again.  When  the  Rayleigh  angle  was  overshot, 
the  arms  were  returned  to  the  nearest  angle  below  it, 
and  the  increment  of  angle  change  was  made  smaller. 

The  overall  behavior  pattern  of  the  minima  was  com- 
plex, so  it  was  frequently  necessary  to  follow  the  be- 
havior of  several  minima  individually  through  a range 
of  angles.  The  result  obtained  would  be  a group  of 
angles , each  of  which  corresponded  to  the  lowest  ampli- 
tude reached  by  one  of  the  minima.  Of  these  angles, 
that  angle  corresponding  to  the  lowest  amplitude  among 
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the  mimimum  values  obtained  was  called  the  Rayleigh 
angle.  In  each  case,  the  lowest  minimum  was  accom- 
panied by  a rapid  variation  of  the  phase  in  its  neighbor- 
hood, amounting  to  as  much  as  180°  of  shift.  This  phase 
variation  provided  some  additional  confidence  in  the 
identification. 

In  early  experiments  the  receiver  and  transmitter 
arms  were  set  at  equal  angles.  Ultimately,  enough 
evidence  was  collected  to  demonstrate  that  the  gonio- 
meter could  not  be  aligned  perfectly  and  that  a genuine 
Rayleigh  phenomenon  could  exist  and  be  identified  with 
the  arms  at  slightly  different  indicated  angles.  As  the 
critical  angle  for  a particular  minimum  was  approached, 
the  arms  could  be  manipulated  independently  and  the 
desired  angle  found  by  watching  the  oscilloscope  CRT 
presentation.  This  technique  produced  angle  settings 
repeatable  within  ±5'  of  arc.  The  difference  in  arm 
angles  at  the  final  settings  are  typically  20'  of  arc,  or 
less.  This  figure  is  consistent  with  estimates  of  the 
alignment  accuracy.  Once  the  arm  angles  were  set, 
amplitude  and  phase  measurement  followed  procedures 
described  in  Secs.  II  and  III.  Amplitude  measurements 
were  made  at  position  intervals  of  0.5  mm. 

V. RESULTS 

The  reflected  intensity  distribution  for  a water— alu- 
minum oxide  interface  is  shown  in  Fig.  5.  The  solid 
curve  is  calculated  from  Eq.  (14).  The  crosses  are 
experimental  points.  The  total  water  path  length  was 
800  mm.  The  incident  beam  width  was  2«)  = 6.35  mm. 
The  incident  angle  a,  was  equal  to  the  Rayleigh  angle 
which  was  14.7°.  The  reflected  field  was  scanned 
across  95  mm,  from  -15  to  +50  mm.  Zero  is  approxi- 
mately the  center  of  the  incident  beam.  The  theory 
predicts  two  minima.  The  second  minimum  was  not 
observable  experimentally  because  the  aluminum  oxide 
specimen  was  not  wide  enough,  but  the  position  of  the 
two  peaks  and  the  first  minimum  agrees  reasonably 
v/ell  with  the  calculation. 

The  stainless  steel  has  smaller  Schoch  displacement 
and  produced  a different  reflected  intensity  distribution 
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FIG.  7.  Reflected  intensity  distribution  for  a water— stainless 
steel  interface  at  *,  = 400  mm. 
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FIG.  8.  Reflected  intensity  distribution  tor  a water-aluminum 
interface. 


as  shown  in  Fig.  6.  The  incident  beam  width  was  6.35 
mm  and  the  total  path  was  800  mm.  The  incident  angle 
was  the  same  as  the  Rayleigh  angle  for  water— stainless 
steel  which  is  30.8°.  A wider  retangular  specimen  of 
stainless  steel  was  used  in  this  experiment.  Now  the 
second  minimum  is  observed  experimentally.  There  is 
also  good  agreement  between  theory  and  experiment. 
The  effect  of  the  changing  of  the  water  path  length 
was  investigated  by  reducing  its  value  to  400  mm.  The 
length  of  each  arm  of  the  goniometer  was  reduced  by 
one-half;  all  other  parameters  remained  the  same.  The 
result  is  shown  in  Fig,  7.  It  is  clear  that  there  is  a 
strong  distance  dependence  on  the  beam  distribution 
which  shows  up  in  the  difference  between  Figs.  6 and 
7.  The  agreement  between  theory  and  experiment  is 
nominally  the  same  at  the  two  distances. 

The  addition  of  z^  dependence  to  the  theory  has  ac- 
counted for  the  diffraction  effect.  In  Fig.  6 the  sec- 
ondary maximum  between  40  and  50  mm  results  from 
diffraction.  Somewhat  less -pronounced  maxima  appear 
in  Fig.  7,  as  well.  Calculation  shows  that  these  sec- 
ondary’ maxima  vanish  as  z^  approaches  zero.  These 
secondary  maxima  also  are  visible  in  schlieren  photo- 
graphs such  as  Fig . 1 . 

Aluminum  was  the  third  material  studied.  Aluminum 
has  the  smallest  Schoch  displacement  of  all  three  ma- 
terials, which  is  16  mm  for  2-MHz  ultrasonic  waves. 
The  Rayleigh  angle  is  30.6°  for  a water-aluminum 
interface.  Both  theory  and  experiment  indicate  (see 
Fig  8)  that  the  reflected  beam  has  only  one  major 
peak.  This  is  quite  different  from  aluminum  oxide  and 
stainless  steel . For  aluminum  the  proper  setting  of 
the  incident  beam  angle  required  careful  alignment. 
Slight  deviation  from  the  Rayleigh  angle  produced 
significant  changes  in  the  intensity  distribution.  For 
example,  a change  of  6'  of  arc  produced  the  intensity 
distribution  shown  in  Fig . 9 . The  position  of  the  prin- 
cipal maximum  has  not  changed  much;  however,  a 
secondary  maximum  is  now  visible  at  28  mm  displace- 
ment, with  another  barely  perceptible  at  40  mm. 
Furthermore,  the  energy  distribution  has  changed  near 
the  specular  beam  and  to  the  left  of  zero  displacement. 


For  the  case  shown  in  Fig.  9 we  also  measured  the 
phase  in  the  reflected  beam  and  compared  with  theory. 

As  can  be  seen,  the  measured  phase  is  in  nominal 
agreement  with  theory.  The  phase  variation  indicates 
that  the  center  of  the  reflected  beam  is  out  of  phase 
by  as  much  as  877  with  respect  to  the  end  points.  This 
phase  plot  is  actually  a measure  of  the  beam  profile 
of  the  reradiated  beam,  indicating  that  the  wavefront 
is  far  from  being  planar  as  is  usually  assumed  in  basic 
boundary-value  problems . 

Finally,  as  a test  of  the  sensitivity  of  the  theory  to 
the  incident  angle,  we  kept  the  distance  z =400  mm  and 
increased  the  incident  angle  by  0.5°  increments.  Two 
calculated  curves  for  0.5°  and  1.0°  greater  than  the 
Rayleigh  angle  are  given  in  Fig.  10.  As  can  be  seen, 
the  theory  also  is  sensitive  to  incident  angle  in  this 
range,  as  it  should  be. 

VI.  SUMMARY 

The  behavior  of  a Gaussian  ultrasonic  beam  incident 
on  a liquid-solid  interface  has  been  studied  in  some 
detail  at  the  Rayleigh  angle— the  angle  at  which  surface 
waves  are  excited  on  the  interface.  The  reflected  beam 
is  displaced  in  the  manner  predicted  by  Schoch;  how- 
ever, the  “Schoch  displacement”  in  general  is  too 
large.  Good  agreement  is  obtained  between  the  experi- 
mental results  and  the  theory  of  Bertoni  and  Tamir 
which  assumes  that  the  incident  beam  couples  resonantly 
into  a leaky  surface  wave  at  the  Rayleigh  angle.  The 
energy  is  reradiated  from  this  leaky  surface  wave  and 
interferes  with  specularly  reflected  energy.  For  water— 


FIG.  9.  Reflected  intensity  distribution  and  phase  profile  for 
a water-aluminum  interface  for  an  angle  of  incidence  differing 
from  9jj  by  6'  of  arc. 
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FIG.  10.  Calculated  intensity  distribution  for  a water-alu- 
minum interface  for  two  angles  near 


AljOj,  which  has  a very  large  Schoch  displacement, 
the  agreement  between  theory  and  experiment  is  fair. 

It  is  better  for  water— stainless  steel,  which  has  a 
smaller  Schoch  displacement.  For  water— aluminum, 
which  has  a relatively  small  Schoch  displacement,  the 
agreement  can  be  quite  good;  however,  sensitivity  to 
incident  angle  makes  necessary  very  careful  alignment 
and  measurement. 

Although  the  approximation  £^  = 0 is  accurate  enough 
in  the  analogous  optical  problem,  the  propagation  dis- 
tance z,  must  be  explicitly  included  in  describing 
ultrasonic  wave  reflection  at  the  Rayleigh  angle.  (We 
believe  that  this  would  also  be  true  for  microwaves 
reflected  internally  from  the  surface  of  a paraffin 
prism. )‘“ 

Thus,  to  give  a complete  description  of  the  phenomena 
observed  at  the  Rayleigh  angle,  one  must  account  for 
three  effects;  (1)  specular  reflection,  (2)  displaced  beam 
resulting  from  leaky  waves,  and  (3)  secondary  beams 


resulting  from  diffraction.  These  three  effects  are 
included  in  the  present  treatment. 
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The  reflection  of  an  ultrasonic  wave  at  a plane 
interface  has  been  studied  in  considerable  detail  in 
recent  years.  If  the  incident  bean  is  infinitely  extended, 
one  can  make  a general  definition  of  the  energy  reflection 
coefficient  and  solve  the  problem  in  closed  form.  In  a 
special  case,  such  as  a water-stainless  steel  interface, 
the  reflected  energy  as  a function  of  the  angle  of 
incidence  behaves  as  shown  in  Figure  1.  Two  critical 
angles  are  obtained.  One  for  the  longitudinal  wave  in  the 
stainless  steel,  and  one  for  the  transverse  wave  in  the 
stainless  steel.  At  these  angles  the  theory  predicts  that 
the  energy  is  totally  reflected.  The  reflection  coefficient 


Figure  1.  Energy  Reflection  Coefficient  for. Ultra- 
sonic Waves  Incident  on  a Plane  Interface. 

has  unit  magnitude.  However,  in  attempting  to  verify  the 
theory  one  typically  obtains  data  such  as  those  shown  in 
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Figure  1.  The  experimental  data  follow  the  theoretical 
curve  for  incident  angles  below  the  critical  angle  for 
shear  waves  in  the  stainless  steel,  29*,  but  above  this 
angle  there  is  considerable  discrepancy  between  the  theory 
and  experiment.  The  theory  predicts  unity,  but  the  experi- 
ment shows  a deep  minimum  at  31*.  The  reason  for  this 
discrepancy  lies  in  the  fact  that  the  actual  ultrasonic 

I -m  . ■ - 

be2uii  IS  not  infinitely  wide.  As  Schoch  has  shown,  the 
reflected  beam  is  "displaced"  down  the  interface  and  does 
not  strike  the  receiving  transducer.  Using  the  theory  of 
Schoch,  one  can  calculate  the  magnitude  of  this  displace- 
ment. The  variation  of  the  displacement  with  angle  of 
incidence  is  shown  in  Figure  2.  The  maximum  occurs  at  31", 
in  agreement  with  the  data  in  Figure  1.  The  reason  for 
this  displacement  is  that  the  incident  beam  is  phase 
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Figure  2.  Schoch  Displacement  as  a Function  of  the 
Angle  of  Incidence. 
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matched  with  a surface  wave  which  propagates  down  the 
interface.  The  angle  at  which  this  occurs  is  often 
referred  to  as  the  Rayleigh  angle  since  such  surface 
waves  were  studied  by  Lord  Rayleigh. 

In  Figure  2 it  can  be  seen,  however,  that  the  calculated 
displacement  is  considerably  greater  than  the  value  measured 
because,  in  fact,  the  energy  is  redistributed  when  it  is 
reradiated  by  the  surface  wave. 

Coupling  with  Surface  Waves  Propagating  in  the  -tx  Direction 
The  situation  can  be  described  in  terms  of  the 
coordinates  shown  in  Figure  3.  The  incident  beam  has 

Z 


incidentX^®r 

beamvX 


f 

reflected 

beams 


. W . . . ..X  . 

SOLID 

Figure  3.  Coordinates  Used  in  Solving  the  Reflection 
Problem. 

coordinates  X.  and  Z.;  the  reflected  beam  X and  Z_.  In 
11  r r 

our  experiments  we  used  an  incident  beam  which  has  a 
2 

Gaussian  profile  as  shown.  The  reflected  beam  then  has 
an  intensity  distribution  that  is  not  simply  a Gaussian 
function  displaced  from  the  origin.  It  is  more  compli- 
cated, as  indicated. 

By  using  the  theory  of  Bertoni  and  Tarair^  one  is  able 
to  redefine  the  reflection  coefficient.  Assuming  a Gaussian 
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incident  beam,  one  finds  a different  reflection  coefficient 
for  each  of  the  partial  plane  waves  in  the  Fourier  trans- 
form representation  of  the  reflected  beam.  In  order  to 
obtain  the  energy  distribution  in  the  reflected  wave,  one 
must  make  a p61nt-by-point  integration  over  the  partial 
waves.  The  result  is  that  with  an  incident  Gaussian  ampli- 
tude distribution 


^inc'’''°5 


ejtP(-(x/w^)  + ik^x] 

/»  w cos9, 
o 1 


{11 


the  reflected  wave  has  an  amplitude  distribution  made  up 
of  two  parts:  a specularly  reflected  component  and  a com- 
ponent reradiated  by  the  surface  wave, 

Vrefi(x,2)  = VqU,z)  + V^Cx,z)  . C2) 
The  specular  component  has  the  form 


expt-(Xj./Wj.)''  + ikj^Xj.] 


/%  w cos6. 
r 1 


The  reradiated  component  has  the  form 

/iTw 

Vj(x^Zj^)  = 2V(j<x^,z^)  II exp(ir^)erfc(YlJ 

8 

where 

T = + i ^ (sinSi  - sinflglw^ 


end 
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(3) 
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(5) 


(6) 


Here  erfc  is  the  complementary  error  function  and  4^  is 
the  Schoch  "displacement*  which  is  characteristic  of  the 
interface.  The  total  reflected  field  at  any  point  is 
given  by  the  (vector)  sum  of  the  two  . complex  numbers 
and  Vj.  The  energy  is  proportional  to  the  product  of  this 
sum  and  its  complex  conjugate. 

The  use  of  this  theory  now  makes  possible  a detailed 
comparison  between  theory  and  experiment  such  as  that 
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shown  in  figure  4.  As  can  te  seen,  the  calculated  curve 
for  a water-stainless  steel  interface  agrees  quite  well 
with  experiraenr.  Data  for  other  solid-liquid  combinations 
we  have  tried  are  in  similar  agreement. 


Displacement  (o.n] 


Figure  4.  Energy  Distribution  in  an  Ultrasonic  Beam 
Reflected  at  the  Rayleigh  Angle. 

Coupling  with  Waves  in  the  -X  Direction 

According  to  the  theory  of  Bertoni  and  Tamir  a surface 
wave  which  propagates  in  the  negative  X -direction  can  be 
generated  under  special  circumstances.  A corrugated  inter- 
face produces  such  a wave . We  have  been  able  to  show 
unequivocally  that  such  a phenomenon  exists.  The  energy  in 
the  reflected  beam  is  "displaced"  in  the  negative 
direction.  This  subject  currently  is  under  investigation. 
(Research  supported  by  the  Office  of  Naval  Research.)  ' 
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An  ultrasonic  beam  incident  on  a liquid-solid  interface  with  superimposed  grating  is  shown  to  be 
displaced  in  a backward  direction  at  certain  angles  of  incidence  predicted  by  the  theory  of  T.  Tamir  and  H. 
L.  Bertoni  [J.  Opt.  Soc.  Am.  61,  1397-1413  (1971)1.  The  analogy  between  the  ultrasonic  phenomenon  and 


the  corresponding  optical  phenomenon  is  cited. 

PACS  numbers:  43.35. -l-d.  68.35. -|-q 

An  ultrasonic  wave,  incident  on  a liquid-solid  inter- 
face at  a certain  angle  of  incidence,  can  couple  into 
leaky  surface  waves,  with  the  result  that  the  reflected 
beam  is  displaced  laterally  down  the  interface.  ‘ 
Recently,  we  studied  this  phenomenon  in  some  detail’ 
and  showed  that  the  theory  of  Bertoni  and  Tamir*  agrees 
well  with  experiment.^  The  theory  of  Bertoni  and  Tamir* 
in  fact  is  an  adaptation  to  the  acoustical  problem  of  a 
theory  previously  used  to  described  the  reflection  of 
light  at  a dielectric  Interface.  ’ In  the  optical  theory, 
Tamir  and  Bertoni®  also  point  out  that  if  the  interface 
has  a periodic  structure  superimposed,  it  is  possible 
to  cause  a leaky  wave  to  propagate  in  the  backward 
direction.  In  such  a case  they  predict  a backward  beam 
displacement.  In  view  of  the  close  analogy  between 
light  reflection  at  a dielectric  Interface  and  ultrasonic 
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wave  reflection  at  a liquid-solid  interface,  we  decided 
to  verify  the  correctness  of  this  prediction  for  ultra- 
sonic waves,  even  though  the  corresponding  optical 
experiment  is  yet  to  be  reported. 

The  situation  is  diagrammed  in  Fig.  1.  A beam  of 
width  2w  incident  at  an  angle  d,  couples  to  a leaky  wave 
propagating  in  the  negative  x direction  indicated  by  the 
heavy  arrow.  The  optimum  angle  for  this  to  occur  is 
given*  as 

where  d is  the  period,  / is  the  frequency,  is  the 
propagation  velocity  in  the  liquid,  and  is  the  propaga- 
tion velocity  of  the  leaky  wave.  At  this  angle  the  incident 

Copyright  © 1976  American  Institute  of  Physics 


456 


456 


[33] 


FIG.  1.  Di.igram  of  incident  beam  coupling  to  a backward-di- 
rected leaky  wave  to  produce  backward  displacement  of  re- 
flected beam. 


beam  is  phase  matched  to  a space  harmonic  of  the  leaky 
wave.  This  leads  to  a reflected  beam  laterally  displaced 
in  the  negative  x direction,  as  indicated. 

The  angle  of  incidence  given  by  Eq.  (1)  is  dependent 
upon  frequency  /,  so  that  for  certain  frequencies  the 
right-hand  side  can  become  greater  than  unity.  This 
corresponds  to  an  imaginary  angle  of  incidence;  i.  e. , 
the  phenomenon  does  not  occur  for  any  real  angle  9,. 

We  used  this  fact  to  produce  the  first  unequivocal 
demonstration  of  this  phenomenon  as  follows. 

Schlieren  photography  was  used  to  visualize  ultrasonic 
beams  in  water  reflected  from  a brass  interface  upon 
which  were  ruled  parallel  grooves  of  period  d = 0. 178 
mm  and  depth  / = 0.  025  mm.  Using  1.  49  x 10® 

cm  sec  and  V^  = 2. 015X10®  cm/sec,  one  predicts  9, 

= 41°  for  6-MHz  ultrasonic  waves.  At  2 MHz,  9, 
becomes  imaginary,  so  the  effect  would  not  be  observ- 
able at  2 MHz. 

A 2- MHz  apodized®  quartz  transducer  was  used  to 
generate  either  6-  or  2-MHz  ultrasonic  beams.  Figure 

2 shows  the  result  for  6 MHz.  The  reflected  beam  is 
displaced  in  the  backward  direction  in  the  manner 
indicated  in  Fig.  1,  because  a leaky  wave  is  excited  in 
the  grating  interface.  (The  direction  of  this  leaky  wave 
suggests  that  a grating  coupler  intended  as  a thin-film 
optical  beam  coupler  would  be  more  efficient  if  it  had  a 
configuration  different  from  that  sketched  in  Fig.  13  of 
Ref.  5.) 

Without  making  any  adjustments  except  changing  the 
frequency  to  2 MHz,  we  obtained  the  photograph  shown 
in  Fig.  3.  The  negative  displacement  is  no  longer 
visible  because  the  leaky  wave  is  not  excited  at  2 MHz. 

Measurement  of  the  angle  of  Incidence  In  Figs.  2 and 

3 results  in  9,  = 22.  5°,  which  is  considerably  smaller 
than  predicted  from  Eq.  (1).  The  reason  for  this  dif- 
ference probably  results  from  the  fact  that  the  value 
V'^  = 2. 015x10®  cm/sec  is  the  velocity  of  a Rayleigh 
surface  wave  on  brass,  whereas  the  quantity  appearing 
in  Eq.  (1)  should  be  the  velocity  of  the  negatively  direct- 
ed leaky  wave,  whose  magnitude  at  present  is  unknown. 

If  one  uses  the  measured  angle  of  incidence  in  Eq.  (1), 
one  can  calculate  the  velocity  of  this  leaky  wave  as 


FIG.  2.  Backwared  displacement  of  6-MHz  ultrasonic  beam  at 
a water-brass  grating  interface. 


FIG.  3.  Ultrasonic  beam  (2  MHz)  Incident  on  water-brass 
grating  interface.  Backward  displacement  is  no  longer 
observed. 
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1.47x10®  cm/sec.  This  means  that  the  velocity  of  these 
leaky  waves  differs  considerably  from  the  velocity  of 
the  Rayleigh  surface  wave  on  a plane  interface.  Further 
experiments  will  be  required  to  determine  whether 
there  is  any  fundamental  significance  to  the  fact  that 
this  velocity  is  close  to  the  velocity  of  ultrasonic  waves 
in  water. 

In  conclusion,  it  is  possible  to  demonstrate  in  a direct 
manner  that  uitrasonic  waves  can  couple  into  a back- 
ward-propagating leaky  wave  if  the  reflecting  interface 
has  a superimposed  periodic  structure.  (This  is  an 
indirect  verification  of  the  validity  of  the  prediction  of 
Tamir  and  Bertoni®  that  light  beams  can  undergo  a 
backward  beam  shift. ) This  fundamental  observation  has 
implications  in  optical  beam  coupler  technology  as  well 
as  surface-acoustlc-wave  technology. 

The  authors  are  indebted  to  Professor  H.  L.  Bertoni 
for  Bringing  this  phenomenon  to  their  attention.  They 
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FORWARD-  AND  BACKWARD-DISRLACEMFNT  OF  ULTRASONIC  WAVES  REFLECTED  FROM 
A WATER-SEDIMENT  INTERFACE 


M.  A.  Breazeale**  and  L.  Bj0rn0 

Fluid  Mechanics  Dept.,  Techn.  Univ.  of  Denmark,  2800  Lyngby,  Denmark 


An  ultrasonic  wave  incident  on  the  interface  between 
water  and  sediment  has  been  studied  by  schlieren 
photography.  Evidence  is  given  that  the  incident 
beam  may  couple  to  backward-directed,  as  well  as  for- 
ward-directed leaky  surface  waves,  with  the  result 
that  the  reflected  beam  may  be  displaced  either  in 
the  forward-  or  the  backward-direction  depending  upon 
frequency  and  grain  size. 


INTRODUCTION 

An  ultrasonic  wave,  incident  onto  a liquid-solid  interface  at  a certain 
angle  of  incidence,  can  couple  with  leaky  surface  waves.  This  coupling 
causes  the  reflected  beam  to  be  displaced  laterally  down  the  interface 
[1].  Recently,  this  phenomenon  has  been  studied  in  some  detail  [2-4] 
and  it  has  been  shown  that  the  theory  of  Bertoni  and  Tamir  [5]  agrees 
well  with  experiments.  The  theory  of  Bertoni  and  Tamir  is  in  fact  an 
adaptation  to  the  acoustical  problem  of  a theory  previously  used  to 
describe  the  reflection  of  light  at  a dielectric  interface  [6].  In  the 
optical  theory,  Tamir  and  Bertoni  also  point  out  that  if  the  interface 
has  a periodic  structure  superimposed,  it  is  possible  to  cause  a leaky 
wave  to  propagate  in  the  backward  direction,  giving  rise  to  a backward 
beam  displacement.  Such  a backward  displacement  of  an  ultrasonic  wave 
reflected  from  a water-brass  grating  interface  was  demonstrated 
recently  ( 7 ] . . 

The  theories  involving  coupling  of  ultrasonic  waves  to  leaky  surface 
waves,  and  their  experimental  verificat?ion,  are  considered  only  for 
liquid-solid  interfaces  which  exhibit  a shear  wave  critical  angle.  A 
shear  wave  critical  angle  is  found  if  > Vg  > V,  v;here  and  Vg  are 

respectively  the  velocities  of  the  longitudinal  and  the  shear  wave  in 
the  solid,  and  V is  the  velocity  of  the  (compress! ona 1 ) wave  in  the 
liquid.  For  such  interfaces,  the  reflection  coefficient  as  a function 
of  the  incident  angle  has  the  behavior  shov;n  in  Figure  1,  in  which  the 
critical  angles  for  the  longitudinal  and  the  shear  waves  in  the  solid 
are  labelled  and  . At  the  angle  labelled  0 , leaky  surface 

waves  can  be  excited  [4]  , and  they  lead  to  a displacement  of  the  re- 
flected beam  as  shown  in  Figure  2.  This  schlieren  photograph  of  a 
1 MHz  ultrasonic  beam  incident  onto  a water-glass  interface  shows  two 
reflected  beams,  the  displaced  beam  and  the  specularly  reflected  beam, 
v^ith  a null  strip  between. 

**On  leave  of  absence  from  Physics  Dept.,  Univ.  of  Tenn.,  Knoxville, 
Tenn.  37916,  U.  S.  A. 


Heretofore,  little  attention  has  been  given  to  experimental  and  LlfeO- 
retical  studies  of  interfaces  for  which  > V > Vg.  In  this  case  one 

does  not  observe  a critical  angle  for  the  shear  wave  in  the  solid,  and 
hence  0 does  not  exist . This  is  the  velocity  relationship  characterizing 

a water-sediment  interface,  since,  typically,  = 1742  m/sec,  V = 1483 
m/sec  and  V^.  = 382  m/sec  [8]-  Even  though  theory  has  not  been  devel- 
oped to  describe  reflection  from  this  class  of  interfaces  in  detail,  we 
report  results  of  schlieren  studies  of  the  reflection  of  ultrasonic 
waves  at  a water-sediment  interface  and  interpret  our  results  by  anal- 
ogy with  results  observed  at  a water-solid  interface. 


THEORETICAL  CONSIDERATIONS 

If  one  assumes  an  infinite  plane  wave  incident  onto  a water-solid  in- 
terface, one  can  calculate  the  reflection  coefficient  from  [9]: 

R = cos^  20g  Zg  sin^  203  - Z 

cos^  20g  + Zg  sin^  20g  + Z 

where  the  normal  impedances  are  defined  by 


p.V  P-V 

7 1 7 = ^ „ ^ .p  (2) 

cos  0.  cos  0^  S cos  0„ 

X L S 

and  the  subscripts  i,  L,  and  S refer,  respectively,  to  the  incident 
wave,  the  transmitted  longitudinal  wave,  and  the  transmitted  shear  wave. 
For  tnree  different  water-sediment  interfaces  these  equations  give  the 
reflection  coefficient  as  a function  of  the  incident  angle  shown  in 
Figure  3.  In  these  curves  we  see  that  there  is  a critical  angle  for 
the  longitudinal  wave,  0^^ , at  which  each  curve  goes  to  unity.  There 

is  not  a critical  angle  for  the  shear  wave  in  the  "solid"  because 
V > V . These  curves  should  be  a good  description  of  the  reflection 

O 

coefficient  if  the  incident  wave  were  plane;  i.  e.  if  the  incident  beam 
were  infinitely  wide. 

In  the  experiments  the  incident  beam  has  a finite  width--as  a matter  of 
fact,  the  incident  beam  has  a gaussian  distribution  across  its  finite 
width  in  order  to  get  rid  of  diffraction  side  lobes.  For  this  incident 
beam.  Figure  3 should  be  a reasonable  approximation  to  the  reflection 
coefficient,  but  there  may  be  detail  in  the  actual  experiment  which 
cannot  be  anticipated  from  Figure  3. 

Since  there  is  a close  analogy  betv/een  light  reflection  at  a dielectric 
interface  and  ultrasonic  wave  reflection  at  a liquid-solid  interface, 
the  prediction  of  a backward  displacement  of  a licht  beam  has  in  fact 
been  verified  for  the  ultrasonic  analogue  [ 7 ].  The  optimum  angle  of  inci- 
dence for  an  ultrasonic  wave  of  beam  width  2w  to  experience  a backv;ard 
beam  shift  by  coupling  to  a leaky  wave  propagating  in  the  negative  x- 
direction,  see  Figure  4,  is  given  by 


®i  = ''[f^  ■ 

v/here  d and  f denote  the  period  and  the  frequency,  respectively, 
while  V and  V are  the  ultrasonic  wave  velocity  in  water  and  the 

f\ 


(3) 


[38] 


propagation  velocity  of  the  leaky  wave,  respectively.  At  this  — 
angle  the  incident  beam  is  phase  matched  to  a space  harmonic  of  the 
leaky  wave,  which- leads  to  a reflected  beam  which  is  laterally  dis- 
placed in  the  negative  x-direction.  Since  the  angle  of  incidence  ‘ 
given  by  Equation  3 is  dependent  upon  both  frequency  and  period, 
and  since  sin  0j^  must  have  values  between  zero  and  unity  for  real 
Gi,  there  is  a relatively  narrow  range  of  grain  sizes  (which  deter- 
mine the  period  d)  for  which  any  specific  frequency  would  be  expec- 
ted to  couple  to  the  leaky  wave. 


EXPERIMENTS  AND  DISCUSSION 

By  using  schlieren  photography  we  have  observed  that  the  reflection 
of  an  ultrasonic  wave  from  a water-sediment  interface  nominally 
behaves  in  the  manner  shown  in  Figure  3 for  0^  < G^l • As  the  angle 
of  incidence  increases  the  reflection  coefficient  increases,  and 
becomes  fairly  large  for  angles  near  to,  and  greater  than,  G^^j^  the 
critical  angle  for  the  longitudinal  wave  in  the  sediment.  Nominal- 
ly, then,  the  reflection  can  be  described  in  analogy  with  a water- 
solid  model,  as  expected. 

However,  even  with  the  qualitative  nature  of  schlieren  observa- 
tions, one  can  see  that  between  G^l  ^nd  90°  the  reflection  coeffi- 
cient does  not  remain  near  unity  as  predicted  by  theory  (see  Fig. 3). 

In  fact,  the  schlieren  pictures  indicate  that  the  reflection  coef- 
ficient in  this  range  of  angles  remains  quite  low  - possibly  one 
tenth  that  predicted  by  theory.  There  are  at  least  two  reasons  for 
the  difference  between  theory  and  experiment.  First,  the  interface 
is  not  sharply  defined  in  the  experiment,  so  that  there  is  a finite 
region  over  which  the  physical  properties  gradually  change  from 
those  of  the  liquid  to  those  of  the  sediment.  Second,  the  sediment 
grains  can  act  scattering  centers  for  the  incident  ultrasonic  wave. 

This  scattering  of  the  incident  wave,  in  addition  to  reducing  the 
reflection  coefficient,  also  can  lead  to  a small  effect  similar  to 
the  backward  displacement  described  above.  If  one  uses  the  mean 
grain  size  of  0.38  mm  in  Equation  3,  and  assumes  that  the  leaky 
surface  wave  excited  at  the  water-sediment  interface  is  of  the  Stone- 
ly  type  with  V|^  = 0.9V,  then  one  can  calculate  that  the  angle  for 
optimum  coupling  with  the  leaky  wave  would  be  Gj^  = 59°  for  2 MHz.  At 
an  angle  of  incidence  of  approximately  57°  one  sees  the  reflection 
phenomenon  shown  in  Figure  5,  in  which  the  striations  in  the  reflec- 
ted beam  presumably  result  from  the  interaction  between  the  specularly 
reflected  beam  and  the  backward-displaced  beam.  These  striations  vanish 
for  larger  or  smaller  angles  of  incidence.  Evidently  the  coupling  to 
the  backward-directed  leaky  wave  is  so  weak  that  a clearly  defined  dis- 
placement as  that  given  in  reference  7 is  not  present  with  this  parti- 
cular water-sediment  interface.  However,  the  photograph  given  in  refe- 
rence 7 does  show  striations  similar  to  those  in  Figure  5. 

Finally,  in  the  range  of  angles  greater  than  the  critical  angle  for 
the  longitudinal  waves  0^^.  , one  observes  reflection  such  as  shown  in 
Figure  6,  in  which  the  reflected  beam  clearly  is  made  up  of  two  parts: 
the  specularly  reflected  beam  and  a displaced  beam,  with  a null  strip 
between  the  two.  Such  phenomena  have  been  studied  in  some  detail  for 
liquid-solid  interfaces  characterized  by  Vj^  > Vg  > V [4].  However, 
the  water-sediment  interface  used  here  is  characterized  by  > V > Vg, 
for  which  presently  available  theory  is  not  applicable.  As  for  as  we 
are  aware,  this  photograph,  taken  at  1 MHz,  is  tie  first  direct  demon- 
stration of  the  fact  that  a forv;ard  displaced  beam  can  be  observed 
slightly  above  the  longitudinal  critical  angle.  The' incident  angle 
here  is  62°,  while  the  calculated  critical  angle  is  = 58.3°. 

Coupling  to  the  displaced  beam  in  the  manner  shown  here  was  observable 


from  approximately  50°  to  angles  somewhat  greater  than  62°.  This  is 
interpreted  to  mean  that  the  fact  that  the  interface  is  granular, 
rather  than  plane,  leads  to  a coupling  with  the  leaky  surface  wave 
over  a considerable  range  of  angles  depending  upon  the  grain  size  di- 
stribution. With  water-solid  interfaces  for  which  similar  phenomena 
have  been  observed  near  the  shear  wave  critical  angle  the  excitation  of 
the  displaced  beam  is  a very  sensitive  function  of  the  incident  angle. 
The  phenomenon  comes  and  goes  within  an  angular  span  of  only  one  de- 
gree of  so  [ 4 ] . 


CONCLUSION 

We  have  studied  the  reflection  of  ultrasonic  waves  from  a water-sedi- 
ment interface,  and  have  found  that  the  reflection  phenomena  nominally 
can  be  described  by  analogy  with  the  reflection  of  ultrasonic  waves  at 
a water-solid  interface.  In  addition,  we  observe  phenomena  which  prob- 
ably are  associated  with  the  backward  displacement  of  the  reflected 
beam  by  the  periodicity  superimposed  on  the  interface  by  the  grains. 

A very  clearly  defined  beam  displaced  in  the  forward  direction  is  ob- 
served over  a considerable  range  of  incident  angles  slightly  above 
the  longitudinal  wave  critical  angle.  Presently  available  theory  will 
need  to  be  modified  to  account  for  this  phenomenon. 

We  conclude  that  in  situations  involving  reflection  from  sediments 
(such  as  those  found  in  underwater  acoustics  experiments) , or  from 
other  interfaces  which  may  have  a superimposed  periodicity,  interpre- 
tation of  results  must  include  the  possibility  that  the  incident  beam 
may  couple  to  either  a forward-directed,  or  a backward-directed  leaky 
wave.  In  sediments  the  density  distribution  at  the  interface  is  not 
as  sharply  defined  as  in  solids.  Nevertheless,  we  observe  coupling  to 
a leaky  surface  wave.  The  angle  at  which  this  coupling  occurs  is  fre- 
quency-dependent for  the  backward-directed  leaky  wave  [7]  but  is  in- 
dependent of  the  frequency  for  the  forward-directed  leaky  wave  [4]. 
These  phenomena  are  to  be  expected  even  for  those  interfaces  charac- 
terized by  > V > Vg.  These  considerations  might  be  important  with 
surface  acoustic  waves  devices  as  well  as  underwater  acoustics. 
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Fig.  1 Reflection  coefficient  for  liquid-solid 
interfaces  characterized  by  V > V > V. 


SEDIMENT  1 

? = 2.03  gm/cm 
V|^=  1836  m/sec 

Vj  = 2S0  m/sec 
9f,=  53.8° 


WATER 

6 = 1,00  6m/cm^ 

V = 1482.9  M/SEC 


SEDIMENT  2 

f = 1.98  gm/cm 
Vl=  1742  m/sec 
= 382  m/sec 
9fr  58.3° 


SEDIMENT  3 
0 = 1.91  gm/cm 
V.  = 1711  M/SEC 


S“  503  m/sec 

ln=  60.1° 


INCIDENT  ANGLE  (DEGREES) 


Calculated  reflection  coefficient  for  three 
different  water-sediment  interfaces. 


Fig.  4 Diagram  of  incident  beam  Fig.  5 Reflection  of  an  ultra- 

coupled  to  a backward-  sonic  beam  from  a v/ater- 

directed  leaky  surface  wave.  sediment  interface. 
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Reflection  of  a Gaussian  Ultrasonic  Beam  fror.i  Itetcr-ricxj Rlas 
(Laszlo  Adler) . 
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Reflection  of  a Gaussian  Ultrasonic  Beam  from 
Water-Plexiglass  Interface 

Laszlo  Adler 

In  most  of  our  previous  studies  for  reflection  of  Gaussian  ultrasonic 
beams  from  liquid-solid  interfaces  we  were  concentrating  on  the  region 
at  or  near  the  Rayleigh  angle.  For  solid  materials  where  the  shear 
velocity  is  smaller  than  the  velocity  of  sound  in  the  liquid  no  real 
value  of  the  Rayleigh  angle  and  of  the  leaky  Rayleigh  velocity  exists. 

Such  a material  is  plexiglass  in  water.  Initial  results  for  reflection 
of  a Gaussian  ultrasonic  beam  from  water-plexiglass  interface  were  taken 
by  a continuous  wave  schlieren  system  and  shown  for  several  angles  of 
incidence  on  Fig.  1.  A double  beam  above  the  longitudinal  critical  angle 
(see  Fig.  1C)  is  observed  which  resembles  the  pattern  one  observes  at  the 
Rayleigh  angle  for  materials  having  real  leaky  Rayleigh  velocity.  The 
reflection  coefficient  calculated  for  an  infinite  plane  wave  from  the 
plexiglass-water  interface  (Fig.  2)  does  not  predict  such  behavior.  The 
appearance  of  the  double  beam  above  the  longitudinal  angle  was  also 
observed  for  the  sediment-water  interface  (see  "Forward-  and  Backward- 
Displacement  of  Ultrasonic  Waves  Reflected  from  a Water-Sediment 
Interface,"  M.  A.  Breazeale  and  L.  BjjJrnd,  Paper  No.  6).  The  appearance 
of  the  double  beam  above  the  longitudinal  critical  angle  is  an  indication 
that  the  finite  beam  couples  to  the  interface  and  that  the  surface  wave 
is  leaking  back  to  the  liquid.  No  theoretical  treatment  of  this  problem 
is  available  at  present.  One  needs  to  investigate  the  behavior  of  the 
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reflection  coefficient  for  finite  beams  at  the  longitudinal  critical 
angle  in  order  to  describe  this  phenomenon  mathematically.  It  may  be 
worthwhile  to  mention  that  this  is  much  more  a true  analog  to  the  optical 
Goos-Ihinchen  effect  (observed  at  the  critical  angle  for  light  waves  from 
dielectrics).  The  reflection  phenomenon  at  the  Rayleigh  angle  is  not 
really  a true  analog  of  the  Goos-USnchen  effect. 


Figure  1 


REFLECTION  OF  A GAUSSIAN  ULTRASONIC  BEAM  FROM 
A WATER-PLEXIGLASS  INTERFACE 
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The  Structure  of  Ultrasonic  Leaky 
Waves  and  Their  Interaction  with 
Subsurface  Flaws 


Abstract 

Recent  experiments  and  theoretical  efforts  indicate  that  the 
interaction  of  a bounded  ultrasonic  beam  with  a liquid-solid 
interface  at  the  Rayleigh  angle  produces  a specular  reflection, 
which  carries  energy  hack  into  the  liquid  from  the  point  of 
incidence,  and  a leaky  Rayleigh  wave,  which  propagates  along 
the  interface  and  continuously  radiates  into  the  liquid.  Because 
the  leaky  wave  penetrates  only  a few  wavelengths  below  the 
solid  surface,  its  potential  for  near-surface  flaw  detection  has 
been  investigated.  An  immersible  acoustic  goniometer  system 
was  used  to  measure  the  amplitude  and  phase  distribution 
across  the  reflected  beam  from  a water-metal  interface  for  an 
incident  Gaussian  beam.  This  amplitude  and  phase  distribu- 
tion, which  had  been  successfully  compared  to  a theoretical 
model,  was  found  to  depend  on  the  following  test  parameters: 
frequency,  distance,  beam  width,  and  material  properties. 
Specimens  with  and  without  artificial  defects  were  investi- 
gated. The  reflected  I reradiatedt  field  is  altered  when  defects 
are  present,  as  shown  by  these  goniometer  measurements  and 
by  schlieren  photography. 


INTRODUCTION 

Surface  (or  Rayleigh)  waves  have  been  studied  in  detail  and 
applied  to  naany  NDT  problems.'  They  are  most  often 
launched  and  received  by  one  or  more  contact  transducers. 
The  "leaky”  wave,  so-called  by  Bertoni  and  Tamir,*  has  a 
surface  component  with  some  properties  similar  to  those  of 
Rayleigh  waves.  However,  it  is  the  differences  between  the 
leaky  wave  system  and  ordinary  Rayleigh  waves  which  make 
leaky  waves  of  potential  value  for  NDT  applications. 

Theoretical""'*  and  experimental"’^  investigations  have  es- 
tablished the  stnicture  and  properties  of  the  leaky  wave  sys- 
tem and  provided  methods  for  computing  the  phase  and 
amplitude  of  the  field.  Simple  field  configurations,  such  as  the 
Gaussian  incident  beam,  can  be  solved  by  analytical  approxi- 
mation; more  difficult  problems  may  require  complete  ma- 
chine computation  of  the  solution. 

The  leaky  wave  system,  generated  by  impingement  from  a 
liquid  onto  a liquid-solid  interface,  includes  a specular,  or 
directly  reflected  component,  and  a surface  component  which 
radiates  (leaks)  energy  back  into  the  liquid.  The  behavior  of 
the  surface  component  depends  strongly  on  the  acoustic  prop- 
erties of  the  interface  materials,  and  the  interference  between 
the  components  in  the  liquid  coupling  medium  produces  the 
unique  features  of  the  system. 
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THEORY 

Reflection  at  a Liquid-Solid  Interface 

Lord  Rayleigh*  first  demonstrated  theoretically  that  sur- 
face waves  can  be  launched  along  a solid-gas  or  solid-vacuum 
interface  by  the  impingement  of  waves  from  the  solid  at  a 
characteristic  angle,  later  named  the  "Rayleigh  angle.” 
Schoch*  studied  the  impingement  of  a beam  from  the  liquid 
onto  a liquid-solid  interface;  he  observed  an  apparent  lateral 
displacement  of  the  reflected  beam  at  a characteristic  angle 
which  he  also  called  the  Rayleigh  angle.  Schoch  also  intro- 
duced the  treatment  of  the  bounded  (i.e.,  finite-width)  beam  as 
a superposition  of  infinite  (width)  plane  waves.  His  treatment 
formed  the  basis  for  the  Bertoni-Tamir  theory"  in  which  the 
reflected  field  is  represented  by 

X 

1 v„n(x,z)=y'  V(k,)R(k,)exp[i(k,x-k,z)ldkj. 

. X 

The  coordinate  axes  are  shown  in  Figure  1;  the  k’s  are  wave 
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Figure  I — Coordinate  system  and  schematic  illustration  of  a 
Gaussian  beam  incident  on  a liquid-solid  interface  at  the 
Rayleigh  angle. 


number  components  in  the  coordinate  directions;  Vik^)  is  the 
angular  distribution  of  infinite  plane  waves,  and  Rik^)  is  the 
interface  reflection  coefficient. 

In  their  development  of  an  analytical  approximation  to  the 
reflection  integral  ( 1 ),  Bertoni  and  Tamir  show  that  it  divides 
into  two  parts  and  can  be  written  as 

2 Vn.||(X,Z)=V|,(X,Z)+ V|(X,Z). 

v„  (x,z)  represents  a specular  reflection,  which  resembles  the 
incident  beam  in  its  amplitude  distribution  but  is  shifted  180° 
in  phase.  v,(x,z)  is  the  surface  component,  which  is  in  phase 
with  the  incident  beam  over  part  of  the  interface  and  out  of 
phase  over  the  remainder;  its  amplitude  distribution  is,  in 
general,  significantly  different  from  that  of  the  incident  beam 
and  the  specular  component. 

For  a Gaussian  incident  beam,  as  used  in  the  experiments 
described  below,  Bertoni  and  Tamir^  obtained  an  analytical 
approximation  valid  at  the  interface.  In  order  to  compare 
theory  with  experiment,  Breazeale,  Adler,  and  Scotf  cor- 
rected the  approximation  for  points  in  the  liquid  halfspace. 
The  expressions  for  the  leaky  wave  field  components  are; 

3 Vii(XpZ,)=  1 exp(  -(Xr/Wr)^  + ik[XrSinH„  + (Zr-Z,i)  cos6)|,]| 

27r  VirWrCOS#,, 

4 V|(x,,Zr)  = -2V|(Xr,z,)  Bxp  (y'*)  erfc  (y)] 

As 

where 

5 W,  Xr 


and 

« w^  = w[l  + 2i(Zr-Z,i)  ]>«, 

kw'''cosW|, 

The  beam  halfwidth  w is  measured  at  z,,,  is  the  liquid- 
solid  equivalent  of  the  Rayleigh  angle,  and  As  is  the  so-called 
"Schoch  displacement.”  As  was  derived  by  Schoch  in  his  origi- 
nal lateral  displacement  theory"  and  was  shown“  to  be  math- 
ematically equivalent  to  a surface  wave  decay  constant  occur- 
ring in  the  Bertoni-Tamir  approximation.  It  is  a complex 


[51] 

function  of  the  acoustic  velocities  ani^  densities  of  the  interface 
media. 

I 

Modeling  Subsurface  Flaw  Interactions 

The  surface  component  of  the  leaky  wave  system  decays 
exponentially  below  the  solid  surface  similarly  to  the 
Rayleigh  wave.  Therefore,  it  is  sensitive  to  subsurface  flaws  in 
the  same  region  of  material,  about  one  Rayleigh  wavelength 
below  the  solid  surface.' 

The  interaction  with  a flaw  was  modeled  by  writing  the 
reflected  field  as 

7 Vr,.||(X,Z)  = V,i(X,z)-l-a(X)V|(X,Z). 

The  defect  tested  was  a side-drilled  hole,  represented  as  a 
cylinder  of  radius  R,  located  at  x, . Then  one  can  write 

1,  for  x<x.  -R 

e a(x)=|  (x-x, -t-R),  for  x,  -R£x£x, +R 

I 2R 
a,  for  X ^x,  + R 

where  a is  the  fraction  of  surface  component  amplitude  re- 
maining after  interaction  with  the  flaw. 

EXPERIMENTAL  APPARATUS  AND 
PROCEDURES 

Transducer  Design  and  Operation 

For  the  present  experiments  it  was  necessary  to  have  a 
well-defined  beam,  free  of  side  lobes.  The  effective  width  was 
not  critical,  although  the  theory'  predicts  some  potential  dif- 
ficulties with  very  narrow  beams.  Gaussian  amplitude  varia- 
tion across  the  wavefront  was  desired  because  this  case  had 
been  solved  analytically,  providing  a good  opportunity  for 
careful  comparison  of  theory  and  experiment. 

For  the  investigation,  the  transducers  (2)  were  fitted  with 
1 in,  diameter  x-cut  quartz  crystals,  plated  with  conducting 
material  (gold  or  silver)  on  the  liquid  side.  The  rear  electrode 
was  a metal  strip  aligned  to  be  symmetric  about  a diameter  of 
the  crystal,  parallel  to  the  goniometer  axis.  The  transmitter 
was  driven  by  a burst  of  2 MHz  RF,  about  6 g-sec  long,  which 
produced  a pulse  in  the  receiver  about  20  gsec  long. 

Results  showed  that  a 3.68  mm  receiver  electrode  consis- 
tently produced  the  narrowest  beam  profiles,  so  that  particu- 
lar size  was  used  in  the  receiver  for  all  experiments. 

A 6.35  mm  electrode  was  used  in  the  transmitter  because  its 
beam  was  narrower  and  diverging  less  rapidly  at  the  distances 
selected  for  reflection  measurements,  and  because  it  produced 
larger  signal  amplitudes  at  the  receiver.  Figure  2 is  the 
amplitude  and  phase  profile  made  with  these  electrodes. 

The  Acoustic  Goniometer 

Figure  3 shows  the  goniometer  removed  from  its  water 
tank.  For  operation  it  is  immersed  in  water  up  to  a point  about 
halfway  between  the  top  ofthe  sped  men  holder  and  the  base  of 
the  column  on  which  the  arms  mount.  The  left  arm  carries  the 
transmitter  and  the  right  arm  carries  the  receiver  (the  slide 
moves  the  receiver  perpendicular  to  the  arm).  The  transducer 
carriages  can  be  moved  radially  along  the  arms  and  clamped 
at  radii  of  about  5 to  50  cm. 

The  specimen  holder  accommodates  rectangular  blocks  (2 
by  4 in.)  up  to  1 in.  thick  and  right  circular  cylinders  (3  in.  in 
diameter)  up  to  a few  inches  long.  Either  specimen  type  is  re- 
strained so  that  its  front  (reflecting)  face  coincides  with  the 
rotation  axis  of  the  arms.  Figure  3 shows  a cylindrical 
aluminum  oxide  specimen  in  place. 

Acoustical  Measurements 

Amplitudes  were  measured  by  reading  the  voltage  devel- 
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Figure  2 — Acoustic  goniometer. 


Figure  3 — Amplitude  and  phase  of  the  acoustic  beam  profde  of 
a Gaussian  transducer. 


oped  across  the  receiver  crystals  from  an  oscilloscope  cathode 
ray  tube. 

Several  potential  applications  motivated  efforts  to  measure 
phase:''  (1)  goniometer  alignment;  (2)  tracing  wavefront  con- 
tours in  transducer-beam  profiles;  and  (3)  tracing  wavefront 
contours  in  Rayleigh  angle  reflections.  A technique  devised 
for  detecting  a constant  phase  led  to  improvements  which 
allowed  approximate  but  repeatable  measurements  of  rela- 
tive phase.  Triggering  the  oscilloscope  from  the  transmitter 
pulse  allows  relative  phase  measurements. 


-“r 1 1 I 1 1 I 
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Figure  4 — Intensity  distribution  of  a reflected  ultrasonic  beam 
from  a water-stainless  steel  interface  at  400  mm  water  path. 


In  this  condition,  the  displayed  signal  from  the  receiver  is 
locked  in  phase  with  the  transmitter  driving  signal.  All  of  the 
pulses  are  coherent  with  one  another  because  they  originate 
from  the  same  CW  source,  and  the  oscilloscope  trigger  circuits 
fix  the  phase  point  on  the  RF  waveform  at  which  the  sweep 
starts.  The  only  remaining  phase  changes  result  from  changes 
in  the  propagation  delay  of  the  acoustic  signal  along  its  water 
path,  and  these  are  observable  on  the  CRT  as  movement  of  the 
entire  waveform  across  the  screen. 

Rayleigh  Angle  Identification 

Reports  of  other  experimenters  established  a range  of  val- 
ues where  the  Rayleigh  angle  for  a material  was  expected  to 
lie.  Starting  at  angles  1-2°  below  the  expected  range,  the 
transmitter  and  receiver  arms  were  set  at  equal  angles.  Then 
the  receiver  was  scanned  through  the  full  range  of  the  cross 
slide,  and  the  locations  and  values  of  all  local  maxima  and 
minima  in  the  amplitude  were  noted.  The  slide  was  set  to  a 
position  corresponding  approximately  to  the  center  of  the  in- 
cident beam.  The  oscilloscope  was  set  to  indicate  zero  relative 
phase  at  this  position,  and  relative  phase  was  observed  as  the 
slide  was  scanned  to  the  right  and  left.  Locations  of  large  or 
sudden  changes  in  phase  were  noted. 

As  the  Rayleigh  angle  was  approached,  the  maxima  and 
minima  occasionally  shifted  position  and  relative  value 
slightly,  and  one  minimum  would  become  an  absolute 
minimum;  i.e. , its  amplitude  would  decline  to  a value  less  than 
those  of  the  other  minima. 

The  overall  behavior  pattern  of  the  minima  was  complex,  so 
it  was  frequently  necessary  to  follow  the  behavior  of  several 
minima  individually  through  a range  of  angles.  That  angle 
corresponding  to  the  lowest  amplitude  among  the  minimum 
values  obtained  was  called  the  Rayleigh  angle.  In  each  case, 
the  lowest  minimum  was  accompanied  by  a rapid  variation  of 
the  phase  in  its  neighborhood,  amounting  to  as  much  as  1 80°  of 
shift.  This  phase  variation  provided  some  additional  confi- 
dence in  the  identification. 

Comparison  of  Reflection  Theory  and  Experimetn 

Amplitude  and  phase  measurements  were  made  on  solid 
stainless  steel  blocks  with  total  water  paths  (transmitter- 
target-receiver)  of  400  and  800  mm  to  observe  the  effect  of 
distance.  The  smooth  curves  compared  with  the  measured 
intensity  data  (computed  by  squaring  the  measured 
amplitudes)  in  Figures  4 and  5 are  computed  values  of 

• |v„.ri|’  = |v„*-V||*. 

The  agreement  between  experiment  and  theory  is  excellent. 
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Figured — Intensity  distribution  of  a reflected  ultrasonic  beam 
from  a water-stainless  steel  interface  at  800  mm.  water  path. 


Figure  6 — Schlieren  photograph  of  the  incident  (left)  and  re- 
flected field  at  the  Rayleigh  angle. 


The  Effect  of  Near-Surface  Flaws  on  Leaky  Waves 

A 6.35  mm  (diameter)  hole  was  side-drilled  0.75  mm  (about 
one-half  shear  wavelength)  beneath  the  2x4  in.  surface  of  a 
stainless  steel  block  to  form  an  artificial  flaw.  To  illustrate  the 
interaction  between  the  leaky  waves  and  the  flaw  a series  of 
schlieren  photographs  was  taken.  In  Figure  6 the  incident 
beam  hits  the  surface  at  a point  far  from  the  flaw  (the  pointer 
indicates  the  positions  of  the  flaw).  The  two  reflected  beams 
correspond  to  the  specular  component  and  the  reradiated 
leaky  wave.  The  null  point  between  these  is  due  to  interfer- 
ence. As  the  incident  beam  is  moved  toward  the  flaw,  the  leaky 
wave  which  penetrates  into  the  .steel  is  interacting  with  the 
flaw  and  is  attenuated.  In  this  case  only  the  specularly  reflec- 
ted beam  shows,  as  in  Figure  7.  When  the  incident  beam  is 
moved  farther  to  the  right,  the  leaky  wave  is  not  attenuated  as 
much  by  the  flaw  and  the  two  beams  again  will  be  visible,  as 
shown  in  Figure  8.  With  this  simple  experiment  we  have  dem- 
onstrated that  near-surface  flaws  can  be  detected  by  leaky 
Rayleigh  waves.  The  experiment  was  repeated  with  a flaw 
drilled  about  two  shear  wavelengths  below  the  surface  and  no 


Figure  7 — Schlieren  photograph  of  the  incident  and  reflected 
field  at  the  Rayleigh  angle.  The  leaky  wave  is  highly  attenuated 
by  the  near-surface  defect. 


Figure  8 — Schlieren  photograph  of  the  incident  and  reflected 
field  at  the  Rayleigh  angle. 


interaction  of  leaky  waves  with  the  flaw  was  observed. 

To  obtain  more  quantitative  information  about  the  varia- 
tion of  the  reflected  field  when  near-surface  flaws  are  present, 
the  goniometer  experiment  was  repeated  at  400  mm  water 
path  and  the  result  is  shown  in  Figure  9.  The  solid  line  is 
calculated  from  Equation  7. 

CONCLUSIONS 

The  agreement  between  theory  and  experiment  for  stain- 
less steel  indicates  that  the  leaky  wave  component  contributes 
strongly  to  the  total  radiated  field  of  ultrasonic  waves  from 
water-metal  interfaces.  The  intensity  distribution  of  the  re- 
flection depends  on  the  ultrasonic  test  parameters  such  as 
water  path  length,  beam  width,  and  frequency  as  well  as  on 
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Figure  9 — Intensity  distribution  of  a reflected  ultrasonic  beam 
from  a stainless  steel  block  with  6.35  mm  hole  located  6.4  mm 
from  incident  beam  center. 


6.  Schoch,  k.Erg.  der  exakten  Naturu  issenschaflen  22  (1950):  160. 

7.  Hildebrand.  B.  P.,  and  Becker,  F.  L.  "Ultra.sonic  Holography  at  the 
Critical  Angle."./  Acoust.  Six-.  Amer.  ,'56  (19741:  4.'59. 


material  properties.  The  leaky  wave  surface  component  which 
penetrates  into  the  solid  is  affected  by  the  presence  of  near- 
surface discontinuities.  It  is  demonstrated  by  schlieren 
method  and  by  goniometer  measurement  that  the  reflected 
field  is  strongly  altered  if  a flaw  is  less  than  approximately 
two  shear  wavelengths  below  the  surface.  This  effect  suggests 
the  possibility  of  using  Rayleigh  reflection  methods  to  detect 
near-surface  defects. 


Application  to  NOT  Situations 

Knowledge  of  the  amplitude  distributions  provides  a basis 
for  the  design  of  near-surface  flaw  inspections.  More  detailed 
studies  of  this  type  would  tell  the  designer  the  optimum  loca- 
tions of  his  transducers  for  maximum  sensitivity  and  allow 
precise  location  of  flaws  with  respect  to  the  transducer  posi- 
tions. The  studies  displayed  here  show  that  neither  the  center 
of  the  specular  reflection  nor  the  classical  "Schoch  distance” 
provides  the  optimum  location  for  the  receiver. 

This  work  also  shows  that  significant  changes  in  phase 
behavior  occur  when  defects  are  present.  The  application  of 
electronic  phase  measurement  to  defect  detection  would  rep- 
resent a significant  step  forward  in  the  use  of  acoustic  infor- 
mation for  materials  inspection. 
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APPLICATION  OF  ULTRASONIC  LEAKY  WAVES  IN  NOE* 

I.aszlo  AJlcr  tinci  CvralH  W.  Scott 
Department  of  Physics,  University  of  Tennessee 
Knoxville,  Tennessee 

ABSTRACT;  Theoretical  investigations  by  Bertoni  and  Tamir  [J.  Appl . Phys,  £,  157  (197J)]  indicate  that  a 
bonndcd~in trasonic  beam  incident  at  the  Rayleigh  angle  to  a linuid-solid  interface  will  produce  a specularly 
reflected  wave  and  a "leaKy  wave"  which  propagates  along  the  interface  and  continuously  radiates  energy  into  the 
liquid.  Since  these  leaky  waves  penetrate  less  than  a few  wavelengths  into  the  solid,  they  should  be  attenuated 
and  scattered  by  the  presence  of  near-surface  defects.  An  immersible  goniometer  system  was  used  to  measure  the 
amplitude  distribution  across  the  reflected  field  from  a water-stainless  steel  interface.  Measurements  indicate 
that  the  presence  of  near-surface  defects  will  alter  the  distribution  of  the  reflected  field.  Similarly, 
alteration  of  the  theoretical  calculation  for  the  reflected  amplitude  is  obtained  when  an  attenuation  mechanism 
is  assumed  for  the  leaky  wave  propagation. 


INTRODUCTION 

The  application  of  surface  or  Rayleigh  waves  to 
NDE  problems  has  been  known  and  well  investigated  for 
some  time.l  Host  often  a contact  transducer  launches 
the  waves  and  receives  echoes  returning  from  cracks  or 
other  defects.  The  technique  described  herein  offers 
some  alternatives  to  conventional  surface  wave 
techniques  by  generating  the  surface  wave  through  fluid 
coupling  and  using  two  transducers  in  a pitch-catch 
arrangement.  It  has  been  recently  established  by  theo- 
retical^ and  experimental^  investigations  that  when  a 
well-defined  (e.g.,  gaussian)  ultrasonic  beam  falls 
onto  a liquid-solid  interface  the  reflected  field  has 
two  components:  (1)  a specular  reflection  and  (2)  a 
surface  or  "leaky  wave"  which  propagates  along  the 
interface  and  radiates  (leaks)  energy  back  into  the 
liquid.  The  total  field,  which  is  the  algebraic  sum  of 
the  two  comjioncnts,  can  be  measured  and  calculated  in 
terms  of  the  ultrasonic  parameters  and  material  proper- 
ties of  the  interface,  Tlie  information  obtained  lends 
itself  to  a quantitative  technique  which  provides 
sensitivity  to  subsurface  flaws. 

THEORY 

Schoch^  was  first  to  study  the  impingement  of  an 
ultrasonic  wave  from  a liquid  to  a solid  interface  and 
observed  an  apparent  lateral  displacement  of  the  beam 
at  an  angle  which  he  referred  tc  as  a Rayleigh  angle. 
Schoch  also  introduced  the  concept  of  treating  a 
bounded  beam  as  a superposition  of  infinite  plane  waves. 
The  reflected  wave  is  represented  by 

\en  ‘ f V{k^)R(k^)exp(i(k^x-k^z)]dk^  (1) 

where  R(k  ) is  the  reflection  coefficient  and  V(k  ) 
depends  on  the  beam  distribution  of  the  incident  beam. 

Schoch  approximates  the  reflection  integral  for  a 
uniform  incident  bean  and  obtains  a formula  which  gives 
lateral  displacement  of  tlic  reflected  beam  but  pre- 
serves an  amplitude  profile  similar  to  that  of  the 
incident  beam.  Experimental  results  of  Breazeale, 
Adler,  and  Smith-  showed,  however,  that  the  reflected 
beam  has  a different  profile  even  when  the  incident  beam 
is  a gaussian  as  illustrated  on  Figure  1. 

Bertoni  and  Tainir^  developed  a new  approximation 
for  the  reflection  integral  (Eqn.  1)  and  derived  a 
solution  for  the  gaussian  incident  beam  which  is  valid 
at  the  interface.  Adler,  Breazeale,  and  Scotf'  au'ded 
correcliotis  to  the  Bertoni -Tamir  solution  for  dist.ance 
from  the  interface.  Their  solution  is  written  as 
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The  beam  halfwidth  w is  measured  at  Zq,  6-  is  the 
liquid-solid  equivalent  of  the  Rayleigh  angle,  Ag  is 
the  so-called  Schoch  displacement,  and  it  depends  on 
the  velocities  and  densities  of  the  two  media.  The 
component  Vq(x  z ) is  the  specular  reflection  whose  pro- 
file rescmbl es''^lhe  incident  wave  and  Vj(x  z^)  is  the 
"leaky  wave"  component  which  penetrates  .about  one  wave- 
length below  tVie  interface.  If  a subsurface  void  or 
flaiv  is  present,  the  leaky  wave  will  interact  with  it. 


Fig.  1.  Schematic  Illustration  of  .a  C.iiissi.in 
Ultrasonic  Bean  Incident  on  a Liquid-Solid  Interface  .>1 
the  Rayleigh  Angle.  Typical  Intensity  laltcrns  of  the 
Reflected  Beams  Arc  Also  Shown. 

MODELING  SUBSURFACE  FlAW  INIF.RACTION 

Let  us  assume  a cylindrical  void  with  radius  R 
whose  center  is  at  x^.  The  reflected  field  is  roditieJ 
as 


and 
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v,cn(='«^)  " * «(x)Vj(..t) 


where 


a(x) 


1 for  x<x^-R 

1 - (x-x^*R)  for  x^-RSxSx^^R 

a for  x>x^*R 


(7) 


(8) 


where  a is  tlic  fraction  of  surface  wave  amplitude 
remaining  after  interaction  with  t)ie  void.  Equations 
(2)  and  (7)  liave  been  experimentally  verified, 

REFbECTED  FIEI.D  MEASUREMENT 


Figure  2 shows  the  goniometer  system  used  to 
measure  the  reflected  field.  The  ttansmitter  is  sus- 
pended from  the  left  arm.  The  receiver  is  suspended 


Fig,  2,  Goniometer  Used  in  Mahing  Intensity 
Measurements, 

from  tlie  cross-slide  mounted  on  the  right  arm.  The 
cross-slide  carries  the  receiver  along  a path  which  is 
perpendicular  to  the  goniometer  arm.  The  transmitter 
was  driven  by  a 2-fniz  pulsed  KF  source.  The  receiver 
output  was  displayed  on  the  delayed  sweep  of  an  oscil- 
loscope for  amplitude  and  relative  phase  neasureraents. 
In  the  experiment  both  transmitter  and  receiver  were 
specially  designed  2-!Clz  gaussian  transducers.  The 
eraplitudc  and  phase  profile  of  the  acoustic  beam  is 
shown  on  Figure  3 as  observed  <t('0  mia  away  in  water. 

The  amplitude  distribution  fits  a gaussian  function 
well.  For  reflection  experiments  both  receiver  and 
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observing  the  RF  display  on  the  oscilloscope  of  the 
reflected  field  amplitude.  At  that  angle  there  is  a 
so-called  "null  point"  observed.  This  is  the  point 
where  the  specularly  reflected  wave  interacts  with  the 
leaky  wave  with  180*  phase  shift  between  the  two  com- 
ponents as  shown  on  Figure  4,  Tlic  samples  used  were  a 
stainless  steel  block  with  polished  surfaces,  one  with- 
out any  subsurface  defects  and  one  with  a near-surface 
defect  simulated  by  a cylindrical  hole  of  3 mm  radius 
drilled  a half  wavelength  below  the  surface. 


■4 


Fig,  4,  Received  Signal  at  the  Rayleigh  Angle, 

The  two  Components  of  the  Reflected  Field  arc  Separated. 

RESULTS 

Tlie  measured  intensity  (computed  by  squaring 
measured  amplitude)  distribution  of  the  reflected  field 
from  a stainless  steel  block  together  with  the  theo- 
retical calculation  is  shown  on  Figure  5,  The  smooth 


Fig.  5.  Reflected  Intensity  nistribution  from  a 
Kater-Stainlcss  Steel  Interface. 


curves  arc  plots  of 


(9) 


Fig.  3.  Amplitude  and  Phase  Profile  of  the 
Incident  Beam. 


The  total  distance  between  transmitter  and  receive)  was 
400  mm.  The  Rayleigh  angle  for  the  interface  is  30.8°. 
The  origin  of  the  coordinates  corresponds  to  the  crniir 
of  the  incident  beam.  Both  experimental  and  thooreli, 
cui-ves  indicate  that  there  is  a small  peak  at  the  left, 
due  to  the  specular  reflection,  and  a larger  peak  (and 
some  diffraction  effects)  at  the  right,  due  to  the 
leaky  wave  component,  with  a well-defined  "null  point" 
where  the  two  components  are  interfering. 

When  the  experiment  is  repeated  under  the  same 
conditions  but  with  a noar-surface  void  in  the  steel 
block,  the  intensity  distribution  is  changed  signifi- 
cantly, as  shown  on  Figure  6.  For  the  flaw,  aVj  was 
used  in  place  of  Vj  in  Eqn.  (9).  11)6  subsurface  void 

is  positioned  7 mt  from  the  center  of  the  incident 
beam.  The  additional  peak  appearing  between  the 
specular  and  leaky  wave  components  is  very  significant. 
Its  center  .aj'pears  to  give  the  location  of  void  x^.  and 
its  width  the  width  of  the  void  2R.  rurthermorc,  the 
leaky  wave  component  is  highly  .attenuated,  as  can  he 
seen  by  comparing  Figure  S and  Figure  6,  The  same 
effect  can  easily  be  demonstrated  qualitatively  by 
observing  schlieren  pbotogr.aphy  as  on  Figure  7.  On 
Figure  7a  the  incident  and  tcflcctcd  field  is  sbowa)  fn 
a stainless  steel  block.  The  position  of  the  Mibsuif.u  ■ 
void  (marked  by  the  pointer)  is  far  from  the  incident 
beam,  llic  reflected  field  is  not  interacting  with  the 
flaw.  Both  specular  component  nnd  the  leaky  wave 


transmitter  arms  are  positioned  to  the  Rayleigh  angle 
with  transducers  equal  dist/uiees  fiom  the  reflector 
(spcciiicn).  The  Rayleigh  angle  can  be  identified  by 
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components  arc  present  with  the  "null  iioint"  hetween 
them.  As  the  incident  beam  is  moved  toward  the  void  at 
one  point  the  leaky  wave  component  of  the  reflected 
field  is  highly  attenuated  due  to  the  interaction  with 
the  near-surface  void  as  shown  on  Figure  7b. 


Fig.  6,  Reflected  Intensity  Distribution  for  a 
Water-Stainless  Steel  Interface  with  a Subsurface  Void, 


Fig.  7.  Schlieren  Photography  of  the  Incident  and 
Reflected  Ultrasonic  Fields,  (a)  The  "beaky  Wave"  is 
Not  Interacting  with  the  Void;  (b)  The  "Leaky  Wave"  is 
Interacting  with  the  Void,  The  position  of  the  Void  is 
Karkcd  by  the  Pointer, 


with  the  experiment.  From  the  reflected  field 
distribution  one  may  obtain  l<oth  position  and  size  of 
a subsurface  flaw.  The  interaction  of  the  leaky  wave 
with  a void  was  also  demonstrated  by  schlieren 
photography, 
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Conclusions 

When  a gaussian  beam  incident  at  the  Rayleigh 
angle  to  a liquid-solid  interface  is  reflected,  the 
reflected  field  distribution  has  two  interacting  compo- 
nents: (1)  specularly  reflected  waves  and  (2)  surface 

"leaky  waves"  radiating  back  to  the  liquid.  The 
reflected  amplitude  distribution  can  be  described  In 
closed  analytic  form  in  tenrs  of  ultrasonic  test 
parameters  and  the  physical  pio;H.rlirs  of  the  liquid 
and  solid,  A iiear-siii  face  void  can  be  modeled  as  an 
attenuation  riethanisia  and  Incorporated  into  the  analytic 
expression.  The  calculation  agrees  reasunably  well 
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ABSTRACT 

The  pjroblem  of  ultrasonic  wave  reflection  at  a Rayleigh  angle  from 
Al^O^  layer-stainless  steel  in  water  was  investigated.  The  so-called 
leaky  Rayleigh  velocity  was  measured  by  using  an  optical  schlieren 
technique  to  identify  the  Rayleigh  angle.  The  magnitude  of  the  measured 
leaky  Rayleigh  velocity  for  t'ue  coated  surface  is  smaller  than  the  corre 
ponding  leaky  Rayleigh  velocity  for  either  stainless  steel  or  for  AI2O2. 


1 


2 


Recent  experimental  investigations  by  Neubaucr  and  Breazeale,  Adler 

2 3 

and  Scott  based  on  the  theory  of  Bertoni  and  Tamir  have  established  that 

a mode  of  elastic  energy  propagation  can  exist  at  a liquid-solid  interface. 

Excitation  of  this  so-called  leaky  Rayleigh  wave  takes  place  when  an 

ultrasonic  beam  of  Gaussian  distribution  is  incident  at  or  near  the  Rayleigh 

angle  to  the  interface.  The  reflected  ultrasonic  beam  contains  a so-called 

null  strip  which  is  the  result  of  the  180°  phase  difference  between  the 

specularly  reflected  beam  and  the  rcradiated  leaky  Rayleigh  wave.  Figure  1 

is  a schlieren  photograph  of  the  incident  and  reflected  ultrasonic  beam 

from  a water-stainless  steel  interface  below,  at,  and  above  the  Rayleigh 

angle.  The  null  strip  appears  when  the  incident  beam  is  at  the  Rayleigh 

angle.  It  has  been  suggested^  tliat  this  null  strip  could  be  used  to 

identify  and  evaluate  the  Rayleigh  velocity  for  the  interface  by  using 


*"R  sine. 


where  Cj^  is  the  leaky  wave  velocity,  0j^  is  the  measured  Rayleigh  angle, 

and  is  the  velocity  of  sound  in  water.  The  leaky  Rayleigh  velocity 

4 

can  be  calculated  exactly  from  theory  by  solving  for  the  roots  of 
secular  equation 


2 2 2 2 ^2 

C,  “C  ^ C 

L J , 


where  one  solution  for  C is  identified  as  C„.  C and  C,  are  the  shear 

K J L 

and  longitudinal  velocity  of  the  solid  and  and  p are  the  density  of 
water  and  tlie  solid,  respectively. 

Tlie  objective  of  this  paper  is  to  report  experimental  observations 


of  the  reflection  of  an  ultrasonic  wave  from  a semiinfinite  solid  when  a 
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thin  layer  of  solid  is  added  to  its  surface.  This  is  an  important 

problem  for  testing  materials  under  coated  surface,  as  well  as  for  some 

surface  wave  device  applications.  The  solid  layer-solid  in  air  is  treated 

by  Famell  and  Adler. ^ According  to  their  analysis,  one  has  to 

distinguish  between  loading  [C^  (substrate)  > (layer)]  and  stiffening 

[C^  (layer)  > C,j,  (substrate)]  of  the  layer  on  the  substrate.  In  the  case 

of  loading  the  Rayleigh  velocity  is  lowered;  in  the  case  of  stiffening 

the  Rayleigh  velocity  will  increase.  No  theoretical  treatment  of  the 

"leaky  Rayleigh  wave"  propagation  for  the  water-solid  layer-solid  exists 

at  the  present  time.  It  is  hoped  that  the  experimental  observations  which 

\ 

are  described  in  the  next  paragraph  will  stimulate  some  investigators. 

Turning  toward  our  experiment,  the  sample  used  is  a stainless  steel 
block  with  dimensions  2.5  cm  x 5 cm  x 10  cm  and  it  is  coated  with  a .24  mm 
thin  aliijninum  oxide  layer  on  the  largest  surface  (shown  on  Figure  2). 

The  thickness  of  the  layer  is  less  than  the  wavelength  corresponding  to 
the  2 Flliz  ultrasonic  beam  used  in  this  experiment.  The  incident  beam 
had  a Gaussian  distribution  and  a half  width  of  4 mm.  Using  a sclilieren 
optic  technique  to  observe  both  incident  and  reflected  beam  from  tlie 
coated  surface,  the  angle  was  adjusted  until  the  null  strip  appeared, 
indicating  the  presence  of  a leaky  Rayleigh  wave.  This  angle  was  recorded 
to  be  35“.  The  corresponding  leaky  Rayleigh  velocity  for  the  coated 
surface  then  is  2.61  x 10^  cm/sec.  The  leaky  Rayleigh  velocities  were 
measured  from  water-stainless  steel  and  from  water-Al203  surface  in  the 
same  way.  The  results  are  given  together  with  calculated  values  (from 
2)  in  Table  1.  The  agreement  between  measured  and  calculated  values 
for  both  cases  are  in  good  agreement.  It  .appears  that  the  Al^O^  layer 
on  steel  lowers  tlic  value  of  the  Leaky  Rayleigh  velocity.  Tlic  same 
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result  was  obtained  on  stainless  steel  when  a zirconium  oxide  layer 
coating  was  used.  These  results  are  the  opposite  of  what  one  would 
expect  for  the  air-solid  layer-solid  case  based  on  the  theory  of  Parnell 
and  Adler. 
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List  of  Figures 

Figure  1.  Reflection  of  a Gaussian  Ultrasonic  Beam  from  a Water- 
Stainless  Steel  Surface  for  Incident  Angles 

a.  25“ 

b.  30.5“  (Rayleigh  angle) 

c.  40“ 

Figure  2.  Stainless  Steel  Block  with  a 0.24  mm  Thick  Al^O^  Layer 
Coating. 


